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Abstract 

In these lectures I review some basic examples of how the concepts of universaUty and scaling can be used to 
study aspects of the chiral and the deconfinement transition, if not in QCD directly but in QCD-like theories. As an 
example for flavor dynamics I discuss a quark-hadron model to describe the phase diagram of two-color QCD with the 
functional renormalization group. Universal aspects of deconfinement are illustrated mainly in the 2 -i- 1 dimensional 
SU(AO gauge theories with second order transition where many exact results from spin models can be exploited. 
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1. Introduction 

Strongly interacting matter fuels the stars and makes 
up almost the entire mass of the luminous universe. 
The underlying theory of quarks and gluons. Quantum 
Chromodynamics (QCD), completely specifies the in- 
teractions. However, these are so complex and non- 
linear that they have yet to be fully understood. In- 
deed, it is these strong interactions which under nor- 
mal conditions confine quarks and gluons into the in- 
terior of hadrons. Understanding the generation of their 
masses, the confinement of quarks and gluons, the dif- 
ferent phases of QCD at extreme temperatures or den- 
sities and the transitions between them are some of the 
great challenges in physics. At temperatures in famiUar 
units close to 2 x lO'^ K ~ 170 MeV, for example, con- 
fined hadronic matter undergoes a transition into an un- 
confined state which is a nearly perfect fluid, the quark- 
gluon plasma. Degenerate neutron matter presumably 
exists in neutron stars at similarly extreme pressures of 
above 1.6 x 10^^ Pa ~ 10 MeV/fin^ Compared to the 
conditions around us, both are of course beyond any 
imagination. The conditions in the finite temperature 
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deconfinement transition are almost 10 orders of mag- 
nitude hotter than the surface of the Earth and still more 
than 5 orders of magnitude hotter than the center of the 
Sun. They bring us back in the evolution of the early 
universe within 10 jus after the Big Bang, and they are 
being recreated in the heavy-ion collision experiments 
at RHIC and LHC [1]. The pressure inside neutron stars 
for comparison corresponds to the weight of about 100 
solar masses pressing on one square meter. Such con- 
ditions at comparatively moderate temperatures require 
heavy-ion collisions at much lower beam energies, cor- 
responding to center-of-mass energies per colUding nu- 
cleon pair in the GeV range instead of TeV at the LHC. 
This energy range was pioneered at the AGS and it will 
be the goal of the HADES and CBM experiments at 
R41R where penetrating probes such as lepton pairs will 
be used to study the ultradense conditions in the early 
stages of the coUisions, in particular [2, 3]. 

Nevertheless, the same characteristic features are be- 
ing discussed in the QCD phase diagram as those, e.g., 
in the phase diagram of water under more normal condi- 
tions of pressure and temperature. Conventionally one 
plots the QCD phase diagram in the plane of tempera- 
ture over baryon chemical potential. An unconvention- 
aUy conventional pressure versus temperature phase di- 
agram for strongly interacting matter is shown in the 
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sketch of Fig. 1. A number of color superconducting 
and other ordered condensed matter phases are expected 
to exist in the high pressure region perhaps comparable 
to the sixteen or so known crystalline phases of water. 
There is the liquid-gas transition to nuclear matter with 
a critical endpoint at a temperature of around 15 MeV, 
and just as in the phase diagram of water there may be a 
second critical endpoint at a higher pressure where the 
conjectured first-order line for chiral symmetry restora- 
tion ends. But also the possibility of an at least approx- 
imate triple point in the QCD phase diagram has been 
discussed [4], which is a point of three-phase coexis- 
tence as in water at the low pressure end of the freezing 
temperature around 610 Pa. 

In order to understand these main characteristic fea- 
tures in the phase digram of strongly interacting mat- 
ter, it has proven to be very useful to deform QCD by 
not only varying the individual quarks' masses but also 
the numbers of their different flavors and colors. Such 
deformations are particularly useful when they lead to 
second-order phase transitions. Then we can apply the 
powerful concepts of universality, scaling and finite- 
size scaling which provide by now standard tools that 
are straightforward applications of the renormalization 
group from statistical physics. The classic example is to 
consider chiral symmetry restoration with only two fla- 
vors of (nearly) massless quarks. In the chiral limit, the 
low temperature phase is characterized by chiral sym- 
metry breaking with the quark condensate {qq)T as the 
order parameter analogous to the spontaneous magne- 
tization in a ferromagnet. Just as the magnetization, 
the condensate melts with temperature due to thermal 
fluctuations. For low temperatures this is very well de- 
scribed by chiral effective field theory [5] based on a 
non-linear realization of chiral symmetry in terms of 
the (would-be) Goldstone bosons which is determined 
entirely by the geometry of the coset G/H or vacuum 
manifold of the symmetry breaking G ^ H. With mass- 
less quarks at two-loop, for example, 

mT-mo[i-^-H^^^ +o(t')]. (1) 

As temperature increases, however, such asymptotic 
expansions necessarily break down eventually. In con- 
trast, the restoration of chiral symmetry at the analogue 
of the Curie temperature in the ferromagnet is very well 
described by a linear sigma model. The global chiral 
symmetry is G = SU(2)rX SU(2)i, corresponding to in- 
dependent flavor rotations of the right and left-handed 
components of the two massless quarks. This symmetry 
is locally the same as SO(4), the group of rotations in 4 
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Figure 1 : Sketch of the QCD phase diagram from [6]. 

Euclidean dimensions. It is broken spontaneously in the 
ordered phase down to the diagonal, vector-like isospin 
symmetry H - S\J{2)v - SO(3). Therefore, scaling and 
universality predict that the singularities of the thermo- 
dynamic observables at this critical point, the endpoint 
of the first-order line in the plane of external field and 
temperature, are all described by the critical exponents 
in the class of the 0(4) Heisenberg ferromagnet in three 
dimensions [7]. In particular, with the analogue of the 
external magnetic field being the explicitly symmetry 
breaking quark mass niq, and t - T/T^ - 1 the reduced 
temperature, this implies that 

{qqh ~ i-tf, for r -> 0" at niq = , and 
{qq)Tc ~ '^y*' for niq ^ Q , where (2) 
/3 = 0.38(1), 5 = 4.82(1). 

Scaling and hyperscaling relations which hold below the 
upper critical dimension d - A then determine all other 
critical exponents, which is referred to as two-exponent 
scaling. More generally, the behavior of the thermody- 
namic observables in the vicinity of the critical point is 
governed by universal scaling functions. For example, 
the magnetic equation of state here is of the form 

{qq)T^m'J'E{tm-'l'P), (3) 

with £(0) = 1 and E{y) ~ {-yf fory ^ -co [8]. 

The study of this 0(4) universality and scaling in the 
quark-meson model with the functional renormalization 
group (FRG) by now also has a long history [9, 10]. 
There is recent evidence, however, that the true scaling 
window might actually be very small, and that it might 
not include the region of physical pion masses in which 
only an apparent scaling may have been observed [11]. 
Finite volume effects were also investigated in the two- 
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triple line, m^i = 



Figure 2: Sketch of a three-dimensional phase diagram with two light 
quarks in the vicinity of a tricritical point at finite baryon chemical 
potential fig = 3fi, by courtesy of the Director of this Winter School. 

flavor quark-meson model and found to be under very 
good control with finite-size scaling [11, 12]. 

When a third, strange quark is included, the situation 
changes depending on its mass. With three massless or 
nearly massless quarks, the transition must be of first 
order because of the axial Ua(1) anomaly [13]. Near its 
physical value, there is nowadays good evidence from 
lattice simulations for O(A^) scaling [14, 15], or perhaps 
at least apparent scaling. If the critical endpoint exists 
in the QCD phase diagram, i.e., for quark masses at the 
physical point, then it seems certainly no less likely to 
exist also when the two light quark masses are sent to 
zero, in which limit it will inevitably turn into a tricrit- 
ical point separating the second-order 0(4) line from a 
triple line at nig = 0, see Fig. 2. With increasing +m^, 
wing lines of critical endpoints would emanate from 
such a tricritical point where three coexisting phases be- 
come indistinguishable. This is a very common view 
but it is seriously being challenged from studies of the 
corresponding critical surface in the three-dimensional 
phase diagram depending on nig, m, and (fi/T)^ through 
lattice simulations at imaginary chemical potential, with 
-(3/7Tf < inlTf < [16]. 

Other important ways to change QCD are to vary the 
number of colors Nc, to change the quarks' color repre- 
sentation or to replace the gauge group altogether. 

For example, the limit of infinitely many colors A^^. 
has inspired many qualitative descriptions of the QCD 
phase diagram [4, 17, 18]. One interesting aspect of this 
limit is that the baryon density becomes an order param- 
eter for Nc — > in particular, also when the number 
of flavors Nf grows along with A^^, i.e. with the ratio 
Nf/Nc held fixed. In the next section we will discuss 
two-color QCD as another example, which perhaps sur- 
prisingly shares some aspects with the large-A^f. limit. 



however with Nc - Nf - 2. Therefore, if the phases 
of many-color QCD with A^^ ~ Nf oa might be re- 
flected in some approximate sense in the real world, it 
seems quite worthwhile to also consider Nc = Nf = 2. 

When changing the quarks' color representation or 
the gauge group as a whole one often exploits the fact 
that a special role in various ways is being played by 
the center of the group, that is the set of all elements 
which commute with the whole group. The center of 
SU(A^) is Zff = {exp(27r/n/A^) : n = 0, 1, . . . A? - 1), 
the set of the A^th roots of unity. The way the center is 
represented in any given irreducible representation de- 
fines its A^-ality A: = 0, 1 , . . . A^ - 1 . If A: is coprime 
with A^, then the center is represented faithfully. This is 
the case in the fundamental representation and all rep- 
resentations with k = I. The other extreme are the real 
representations with A^-ality k = such as the adjoint 
representation in which the center is represented triv- 
ially, i.e. in a fully degenerate way by the unit element. 
Another special case worth mentioning are representa- 
tions with k = A^/2 for even A^ in which the center is 
represented by the real Z2 = (1,-1). These represen- 
tations are called pseudo-real because they are isomor- 
phic to their complex conjugates. One special interest in 
quarks in real or pseudo-real color representations de- 
rives from the fact that the fermion determinant remains 
positive or real, respectively, even at finite chemical po- 
tential. This means that such theories can be simulated 
at finite baryon density without fermion-sign problem 
on the lattice (using pairs of degenerate flavors, if nec- 
essary in the worst case). 

For example, the center of the SU(2) gauge group of 
two-color QCD is Z2 = (1,-1). Its integer representa- 
tions, i.e. the proper rotations of SO(3), have A^-ality 
k = and are real while the half-odd integer representa- 
tions with k = 1 such as the fundamental representation 
are pseudo-real, and this exhausts all possibilities in this 
case. We will have much more to say about the con- 
sequences in the next section. Here it suffices to make 
one general observation, and that is that color singlets in 
product representations can occur only for even powers 
of pseudo-real representations. Therefore, theories with 
quarks in pseudo-real color representations can contain 
only bosonic but no fermionic baryons. Theories with 
fermionic baryons require singlets in odd powers of the 
quarks' representations which singles out the real ones 
if the sign problem is to be avoided. Adjoint fermions in 
SU(A^) gauge theories are the most obvious candidates. 
Other important examples include the compact excep- 
tional Lie Groups G2, F4 and which all have a trivial 
center. The smallest one of these is G2. It has rank 
2 and dimension 14, it is simply connected and it can 
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conveniently accommodate SU(3) as a subgroup. All 
its representations are real, so there is no fermion-sign 
problem even for a single flavor and lattice simulations 
at finite baryon density are possible [19]. Moreover, if 
the G2 gauge symmetry can be broken down to SU(3) 
via the Higgs mechanism in a controlled way, one can 
study how the sign problem gradually reemerges. 

We should eventually be able to completely under- 
stand the phase diagrams of QCD's closest relatives 
without the fermion-sign problem, and this should be 
a worthwhile exercise also for QCD itself. QCD-like 
theories can be used to study universal aspects but also 
as benchmarks for functional methods which work for 
QCD as well but which rely on truncations or model 
assumptions. While two-color QCD with its bosonic 
baryons has been studied extensively in the past, as will 
be reviewed in the next section, studies of the phase di- 
agram of G2 gauge theory are still in their infancy. They 
wiU be an important next step in understanding a sign- 
problem-free variant of QCD with fermionic baryons. 

The other special role that is being played by the cen- 
ter of the gauge group is via the global center synmietry 
and its spontaneous breaking in the deconfinement tran- 
sition of the pure gauge theory. Again, if we change the 
number of colors or the dimension of space from three 
to two, this transition becomes second order and we can 
gain a very precise understanding of it from universality, 
scaling and finite-size scaUng. This wiU be discussed in 
Section 3. 

2. Quark-meson-diquark model of two-color QCD 

Quantum Chromodynamics with two colors (QC2D) 
has been well studied for many years within chiral ef- 
fective field theory and random matrix theory [20-28 J, 
in lattice simulations [29-36], and the Nambu-Jona- 
Lasinio model [37^6]. In this section we focus on 
the (Polyakov-)quark-meson-diquark model for study- 
ing the phase diagram of QC2D with the functional 
renormalization group [47], including fluctuations due 
to collective baryonic excitations. 

The most important differences between two and 
three colors both follow from the property of the SU(2) 
gauge group of QC2D that its representations are either 
real or pseudo-real. This leads to an antiunitary sym- 
metry in the Dirac operator [21]. As one consequence, 
the fermion determinant remains real for non-vanishing 
baryon or quark chemical potential, 0, as it does 
for adjoint quarks in any-color QCD, or in the G2 gauge 
theory with fundamental fermions. Thus, at least for an 
even number of degenerate fundamental quark flavors in 



QC2D there is no fermion-sign problem and the phase 
diagram is amenable to Monte-Carlo simulations. 

Another consequence of the reality or pseudo-reality 
of the quarks' color representation is a PauU-Giirsey 
symmetry which allows to combine quarks and charge- 
conjugated antiquarks into enlarged flavor multiplets. 
As a result, for vanishing chemical potential and quark 
mass, ju = = 0, the usual SU(A?/) X SU(A?/) X U(1)b 
chiral and baryon number symmetries are replaced by 
an extended SU(2iV/) flavor symmetry. A similarly ex- 
tended flavor symmetry is also known from (even num- 
bers of) fermions in 2 H- 1 dimensions with U(l) or 
SU(A^) gauge fields with N > 3. There it is due to 
the fact that there is no physical helicity and that the 4- 
dimensional spinor representation of Dirac fermions is 
reducible in the 3-dimensional space-time. The result- 
ing extended U(2Nf) flavor symmetry is relevant for the 
expected semimetal-insulator transition in QED3 with 
not too many flavors, estimated for Nf < 4, or for 
the description of the electronic excitations around the 
Dirac points of graphene at half-filling which are also 
described hy Nf - 2 such 4-spinors of 3-dimensional 
Dirac fermions in a single layer or Nf = 4 in a double 
layer, for example [48]. With (pseudo-)real gauge fields 
such as SU(2) or G2 in 3 dimensions one obtains an even 
larger extended flavor synunetry [24]. 

As usual, the extended flavor symmetry of our 4- 
dimensional QCD-like theories is (spontaneously) bro- 
ken by a (dynamical) Dirac mass. The breaking patterns 
are somewhat different from the usual, however. In the 
pseudo-real case of the fundamental two-color quarks 
the (non-anomalous) extended SU(2Af/) gets broken 
down to the {2Nf + \)Nf dimensional compact symplec- 
tic group Sp(Nf). In the real cases of adjoint SU(A^) or 
fundamental G2 quarks, for example, the corresponding 
breaking would be S\J(2Nf) Spin(2A?/), flie double 
cover of the proper rotation group SO(2Nf). 

For Nf - 2 the extended SU(4) flavor symmetry 
group and its Sp(2) subgroup combining the isospin and 
baryon number symmetries of QC2D are locally iso- 
morphic to the rotation groups SO(6) and SO(5), re- 
spectively. The coset is given by 5^, the unit sphere 
in six dimensions, and a spontaneously generated Dirac 
mass will lead to five instead of the usual three Gold- 
stone bosons in this case, the three pions plus a scalar 
diquark-antidiquark pair. 

Moreover, for Nc - 2 these color-singlet scalar di- 
quarks play a dual role, as would-be-Goldstone bosons 
and bosonic baryons at the same time. While this 
thus represents the perhaps most important difference as 
compared to the real world, it also makes it much eas- 
ier to investigate the effects of baryonic degrees of free- 
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dom on the phase diagram in functional approaches. In 
that sense the quark-meson-diquark (QMD) model can 
be considered as a first step towards their inclusion in a 
'quark-meson-baryon' model for real QCD with A/^ = 3. 

For the same reason the QMD model of QC2D pro- 
vides a relativistic analogue of the BEC-BCS crossover 
in ultracold fermionic quantum-gases, which has also 
been described successfully with functional renormal- 
ization group methods [49-51]. In contrast to non- 
relativistic models of the BEC-BCS crossover, an in- 
teresting additional constraint thereby arises from the 
Silver Blaze property [52]: When a relativistic chemi- 
cal potential // is coupled to degrees of freedom with a 
mass gap A, as temperature approaches zero, the par- 
tition function and hence thermodynamic observables 
must actually become independent of the chemical po- 
tential as long as yU < A. In general, it is not guar- 
anteed that this constraint is automatically satisfied in 
non-perturbative approaches such as fimctional renor- 
malization group studies which rely on truncations, but 
it is yet another example of valuable extra information 
to devise intelligent truncation schemes [47]. 

One main virtue of the QMD model, however, is to 
explicitly demonstrate the impact of baryonic degrees 
of freedom on the phase diagram by comparing it with 
the corresponding purely mesonic model, as represen- 
tative of typical three-color QCD model calculations. 
For this comparison it is more appropriate to think of 
the vacuum diquark mass as the baryon mass mb rather 
than the pion mass m,,. In QC2D with its extended fla- 
vor symmetry they are the same, but the essential aspect 
of this assignment is that a continuous phase transition 
at zero temperature occurs at a critical quark chemical 
potential jic = ms/Nc. Except for the scale separation 
between m„/2 and m^/S in the real world, this transi- 
tion is then to be compared to the liquid-gas transition 
of nuclear matter in QCD with three colors which is of 
first order, involves the binding energy, and thus occurs 
somewhat below ju = ms/Nc. 

As temperature increases the liquid gas transition 
ends, turning into a crossover with continuously vary- 
ing but nevertheless probably still relatively abruptly in- 
creasing baryon density along some narrow region. This 
rapid increase is generally expected to lead to the strong 
chemical-potential dependence of the chemical freeze- 
out line observed in heavy ion collisions at center-of- 
mass energies below about 10 GeV per nucleon pair, 
the baryonic freeze-out [53, 54]. One might conclude 
that the phase transition line for diquark condensation, 
where a rapidly increasing baryon density develops, 
would be the origin of a corresponding baryonic freeze- 
out Une in two-color QCD, with = Nf = 2 arguably 



not necessarily further from reality than the large Nc 
limits. As in the latter, one might then even identify a 
two-color version of quarkyonic matter [4, 17, 35, 42]. 

Finally, it is worth noting that the model, the func- 
tional renormalization group equations and the tech- 
niques to solve them have a broad scope of applications 
beyond two-color QCD. One example is QCD with 
two light flavors at finite isospin chemical potential, 
which has been studied with the NJL model in mean- 
field plus random phase approximation (RPA) [55, 56]. 
There is a precise equivalence between the correspond- 
ing quark-meson model with isospin chemical poten- 
tial and the quark-meson-diquark model of two-color 
QCD discussed here. Besides changing Nc this merely 
involves reducing the number of would-be-Goldstone 
bosons from five to three again, retaining only one of 
the three degenerate pions for the neutral one, and rein- 
terpreting the diquark-antidiquark pair of QC2D as the 
charged pions of QCD with isospin chemical potential 
[57]. Similar models are also studied in the context of 
color superconductivity [58-60]. The capacity to nu- 
merically solve functional renormalization group equa- 
tions on higher dimensional grids in field space is gen- 
erally useful for competing symmetries, as in a quark- 
meson model study of the axial anomaly with scale de- 
pendent 't Hooft couplings, for example [61]. 

2. 1. Extended flavor symmetry and model Lagrangian 

As all half-odd integer representations of SU(2), its 
fundamental representation is pseudo-real because it is 
isomorphic to its complex conjugate representation with 
the isometry given explicitly by 5 = icrj, = -1. 
Therefore, charge conjugation of the gauge fields in 
QC2D can be undone by the constant SU(2) gauge trans- 
formation S = i(T2. With T" = cr"/2foT the fundamen- 
tal color generators in SU(2), one has 

T"'^ = 7"* = -ST''S~K (4) 

We will reserve cr, (r,) for the Pauli matrices in spinor 
(flavor) space from now on. Together with the charge 
conjugation matrix C in spinor space, hkewise with 

= -1, and complex conjugation denoted by K one 
then defines an antiunitary symmetry T = SCK with 

= +1 as time-reversal invariance in quantum me- 
chanics which leaves the Dirac operator invariant. 

For comparison, the irreducible representations of the 
proper rotations, as the adjoint representations of all 
SU{N), are examples of real representations with her- 
mitian generators that satisfy 
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More generally this corresponds to an isometry S for 
complex conjugation as in (4) but now with 5^ - +1. 
For all real color representations, as those of the ad- 
joint groups SU(A^)/Ziv or the exceptional Lie groups 
mentioned in the introduction such as G2, one thus has 
= -1, correspondingly. The Dirac operator then 
has an antiunitary symplectic synnmetry which results 
in a two-fold degeneracy of its eigenstates and a positive 
fermion determinant even for a single flavor. This leads 
to the classification of the Dirac operator by the Dyson 
index j3 of random matrix theory [20, 21], with jS - \ 
for fermions in pseudo-real color representations corre- 
sponding to the Gaussian orthogonal ensemble, j8 = 4 in 
real color representations corresponding to the Gaussian 
symplectic, and /3 - 2 corresponding to the Gaussian 
unitary ensemble otherwise. The Dirac operator in the 
fundamental color representation of QC2D (and in the 
A^-ality N/2 representations of the SU(A'')'s with even 
A^) falls in the first class (jS = 1). Those of G2, F4 and 
Eg (and the A^-ality zero representations of SU{N)) fall 
in the second (/? = 4).^ 

Following [21], let's start from the standard kinetic 
part of the EucUdean QC2D Lagrangian, in the chiral 
basis, 

-Ckin = tf>ll>>/f = i/flio-nDutfTL - i//li(rlDuifrR , (6) 

with hermitian y-matrices, cr^ = {-i, &), and iItril, >lfg/i 
as independent Grassmann variables with 1^^^^ = 
The covariant derivative is D^-d^ + iA^, and the cou- 
pling is absorbed in the gauge fields = A^T". 

The two terms in (6) get interchanged under the anti- 
unitary symmetry T. One can apply a corresponding 
transformation to only one of the two terms to change its 
sign, however Using (-ia-2) for the chiral /^-component 
of the charge conjugation matrix C in the second term, 
say, by changing variables to {//r = -icr2Stl/*jf and 1^^ = 
-icr2Si/fR, one can therefore re-express 

Xkin = "V^a^D^V (7) 

in terms of 2Nf 4-dimensional spinors *P = (il/i, ij/Rf 
and = 1^^). Because Xidn is now block diago- 
nal, the SU(2A^/) symmetry in the space combining fla- 
vor and transformed chiral components is manifest in 
this form. With the same transformation of variables 



Without the extra minus sign from the charge conjugation matrix 
C after spin diagonalization, the antiunitary symmetries of staggered 
fermion operators for SU(2) are opposite to those of the continuum 
or other lattice Dirac operators, /J = 4 for fundamental (pseudo-real, 
half-odd integer) and j8 = 1 for adjoint (real, integer colored) quarks. 



the quarks' Dirac -mass term becomes 

miA<A = ^(*P^icr25So>P - >P*^j(r252o*l'*) , (8) 
where the symplectic matrix 

acts in the 2A^/-dimensional extended flavor space and 
transforms as Eq U^'LqU, because it is antisymmet- 
ric according to the Nf(2Nf - l)-dimensional antisym- 
metric rank 2 tensor representation of SUi2Nf). The in- 
variance group of Eq as bihnear form on complex 2Nf- 
vectors is the compact symplectic group Sp(A'/) which 
is the intersection of the unitary \J{2Nf) and the sym- 
plectic Sp(2Nf, C), therefore sometimes also referred to 
as USp(2A'/). An explicit(dynamical) Dirac mass thus 
expUcitly(spontaneously) breaks the original SU(2A^/) 
down to Sp(A^/). 

For real color representations (J3 = 4), going through 
the same steps with replacing 5 — » 1 [21], the mass 
term is a symmetric color singlet and the correspond- 
ing flavor matrix is therefore symmetric, likewise. It 
then belongs to the Nf(2Nf + l)-dimensional symmetric 
rank 2 tensor representation of SU(2A^/), and the invari- 
ance group, generated by the Nf{2Nf-\) antisymmetric 
hermitian 2N f x 2Nf matrices, is SO{2Nf) or its double 
cover Spin(2A^y ) for fermionic states. 

For Nf - 2 flavors the enlarged flavor symmetry 
group is SU(4), not U (4) because of the axial anomaly, 
it replaces the usual chiral and baryon number symme- 
tries SU(2)£ X SU(2)jj X U(1)b. Just as fliis extended fla- 
vor SU(4) shares its 15 dimensional Lie algebra with the 
group of rotations in 6 dimensions, SO(6), its Sp(2) sub- 
group leaving the (/3 = 1) Dirac -mass term invariant has 
the 10 dimensional Lie algebra of SO(5) (in fact they 
are both the universal covers of the respective rotation 
groups). So the symmetry breaking patterns by Dirac 
mass terms are locally the same as SO(6) SO(5) and 
S0(6) ^ S0(4) foryS = 1 and yS = 4, respectively. 

On the other hand, at finite chemical potential //, but 
still with massless quarks, m = 0, the SU(2A^/) symme- 
try is broken explicitly by [21] 

liiftyoilf = juY^Bo^P , with Bo = (^J j , (10) 

down to SU(A^/)l x S\i{Nf)R x U(1)b. For Nj- = 2, 
in terms of the rotation groups, this symmetry breaking 
pattern is locally flie same as SO(6) SO(4) x SO(2). 

When both and m are non-zero, the unbroken flavor 
symmetry is of course given by the common subgroup 
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S\](2)v X U(1)b of the two limiting cases: (i) /i — > 
at finite m with either Sp(2) (J3 - 1) or Spin(4) (fi - 4) 
symmetry, and (ii) m — > at finite fj. with SU(A^/)l X 
SU(Nf)R X U(l)s in either case, as discussed above. 

Whether the resulting SU(2)y x U(1)b is actually 
closer to the combined isospin plus baryon-number 
symmetries or the standard chiral symmetry, naturally 
depends on the relative sizes of the quarks' Dirac mass 
m and their chemical potential fi. More precisely, for 
quark chemical potential jj < m„/2 the chiral sym- 
metry breaking pattern essentially remains the ju = 
one and the vacuum alignment is said to be (^^)-like. 
There is only an explicit breaking of the combined 
isospin/baryon-number symmetry proportional to fi: ac- 
cording to Sp(2) — > SU(2)yxU(l)B which is the same as 
SO(5) ^ SO(3) X SO(2) for the bosonic states of QC2D 
with its pseudo-real fundamental quarks (J3 - 1), or ac- 
cording to Spin(4) = SU(2) x SU(2) ^ SU(2)v x U(1)b 
for G2 with real fundamental quarks or other theories 
with quarks in real color representations (J3 = 4). If 
the chemical potential is small, the respective enlarged 
isospin/baryon-number symmetries are approximately 
realized. The vacuum is unchanged and the degenera- 
cies in the spectrum simply split up by quark num- 
ber proportional to fj.. For example, among the would- 
be-Goldstone bosons the masses ota* of diquark and 
antidiquark in QC2D split from those of the pion as 

When Ota- = 0, i.e., for yu > OT;r/2, a diquark con- 
densate develops at zero temperature and the vacuum 
alignment starts rotating from being (^^)-like to be- 
coming more and more (^^)-like as /i is further in- 
creased. The chiral condensate then rapidly decreases 
and chiral symmetry gets restored to the approximate 
SU(2)z, X SU(2)r ^ SO(4). In this phase, the baryon 
number U(1)b is spontaneously broken and the remain- 
ing isospin SU(2)v changes from the approximately re- 
alized enlarged isospin/baryon-number Sp(2) symmetry 
to becoming an approximate standard chiral S\J{2)i x 
SU(2)r symmetry.^ Both are only approximate sym- 
metries, the first is weakly broken by the small jU (but 
there is a large dynamical quark mass), while the sec- 
ond is weakly broken by the small current quark mass 
OT (at large //). They change in a crossover In fact, 
the dynamical contribution to the quark mass changes 



^In our example with Nf = 2 quarks in a real color representation 
it might seem that there is no change in the approximate symmetry, 
it is Spin(4) = SU(2) X SU(2) for small chemical potential below the 
crossover and SU(2)£ X SU(2)k above. They may look the same but 
they are two different subgroups of the original SU(4), the first one 
contains U(1)b as a subgroup and the second one does not. 



r 
> k 




Figure 3: Schematic phase diagram for QC2D in the parameter space 
of temperature T, quark mass and baryon chemical potential /ig. 



in this crossover from being predominantly the original 
Dirac mass, with a condensate of tightly bound light di- 
quarks, to becoming a dynamical Majorana mass as in 
BCS theory. This vacuum realignment in the diquark- 
condensation phase with superfluidity of the bosonic 
baryons in QC2D is the analogue of the BEC-BCS 
crossover in ultracold fermionic quantum gases. 

At zero temperature, for 2// < m„, below the onset 
of the condensation of diquarks as bosonic baryons, the 
baryon density remains zero and the thermodynamic ob- 
servables must be independent of fi. Because this is far 
from obvious to verify explicitly in actual calculations, 
it has been named the Silver Blaze Problem [52]. In or- 
der to be able to excite any states at zero temperature, 
and with a gap in the spectrum, the relativistic chemi- 
cal potential needs to be increased beyond the mass gap 
in the correlations to which it couples. Here, with a 
continuous zero-temperature quantum phase transition 
at fig = 2fi = m„ this gap is simply given by the lightest 
baryon mass in vacuum which because of the extended 
flavor symmetry in QC2D coincides with the pion mass, 
mB - inn- This latter property is of course special to 
Nc - 2 OT other theories with quarks in pseudo-real 
(without fermionic baryons) or real color representa- 
tions (with fermionic baryons, but where the lightest 
baryon with mass ot^ = m„ is still bosonic). The Sil- 
ver Blaze property must hold as it does here, however, 
also when there are no bosonic baryons as in QCD up 
to a quark chemical potential of the order of nig/Nc (re- 
duced by 1 /Nc of the binding energy per nucleon when 
the transition is of first order). 

At finite temperature, a qualitative picture emerges 
for the phase diagram of QC2D as sketched in Fig. 3. 
The solid line in the T - plane represents the contin- 
uous zero-temperature transition with diquark conden- 
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sation which is of mean-field type. Because the quark 

mass niq scales quadratically with the pion mass, it will 
occur along a parabola niq oc /d^. The thick dashed lines 
represent the corresponding second-order transitions at 
finite temperature in fixed nig planes of the 0(2) univer- 
sality. The thick fine along the temperature axis is the 
magnetic first-order transition in the fiB = plane which 
might end in a multicritical point. When viewed in the 
I^B - plane, this is the critical endpoint in the 0(6) 
universality class for the chiral phase transition in QC2D 
with its extended SU(4) flavor symmetry. In the nig = 
plane, the vacuum alignment will always be (<j'<j')-like, 
for no-matter-how- small fi > 0. Therefore, in this plane 
one only has the second-order 0(2) line which, if it ends 
in the same point, would make it multicritical. 

The construction of the QMD model for QC2D starts 
from the flavor structure of the standard chiral con- 
densate and the quark mass term which is of the form 
*P^Eo*P. It therefore transforms under the full flavor 

SU(4) according to the six-dimensional antisymmetric 
representation in the decomposition 4i8>4 = 6©10. 

The other components belonging to the same multi- 
plet are obtained from transformations 

^'^U'¥ , U = expiiff'X'') € SU(4)/Sp(2) . (11) 

Then, "P^Eq*!' ^^S^P, where, from Cartan's immer- 
sion theorem, the whole coset SU(4)/Sp(2) = is 
obtained in this way via E = U^'LqU. The coset ele- 
ments S are in turn parameterized by six-dimensional 
unit vectors ft d&Y. = nt, with SjSy + Sjs, - 26ij and 

E = (Eo, moX") such that X", a = 1 ... 5, form a ba- 
sis for the coset generators [25]. Thus, one verifies ex- 
plicitly that the vector ^'^I,^ transforms as a (complex) 
six-dimensional vector under SO(6). 

The QMD model for QC2D is therefore defined 
by couphng the real SO(6) vector of quark bilinears 
(-vpr^vp + h.c.) to a vector of mesonic fields <p - 
(cr, TT, Re A, Im A)^ formed by the scalar cr meson, the 
pseudoscalar pions 7? and the scalar diquark-antidiquark 
pair A. This yields the Lagrangian, 

, (12) 
- Y*^j(r25f»I'*) ■$+ -idjf + Vi$), 

where V0) is the meson and diquark potential. A 
non-vanishing chemical potential not only couples to 
the quarks but also to the bosonic diquarks. Rewriting 



Eq. (12) in terms of the original variables, one obtains 

+ I (A*(^^Cr'T25^) + A(il^^'Cy^T2Sr)) 

+ \ i(d, - 2fi ^)A)(5^ + 2fi 6^)A* , (13) 

with C = y^^ and a flavor- and color-blind Yukawa 
coupling g. With 

y(0)= ^(^-vY-ccr, (14) 

one obtains the corresponding 0(6) linear sigma model; 
and in the limit A ^ 00, the bosonic part of -Cqmd is 
then equivalent to the leading-order ;ifPT Lagrangian of 
Refs. [21] with the identifications v - f,t - 2F and 
c = f^ml - 2Fm^. The coefficient of the leading term 
in n of the;(fPT Lagrangian, which is i^ir(LB^Y!^ B) with 
B - UB()U \ was fixed from gauging the flavor SU(4) in 
[20]. Here it simply follows from -2/i^|Ap as part of the 
kinetic term of the complex scalar diquark field A with 
chemical potential hb - 2//. This implies in particular, 
that the meson/diquark potential V{$) itself, up to the 
explicit breaking by -ccr, which in principle needs to 
be only S0(4) x S0(2) invariant at finite fi, must remain 
SO(6) invariant, however, at this leading order, OQj^), 
and therefore at 0(0^) in the fields, Ukewise. We can 
thus only have an SO(6) invariant mass term in V($). 

Gauge field dynamics and confinement effects can 
be modeled also in QC2D by including a constant 
Polyakov-loop variable as a background field as in the 
NJL model [42], and analogous to what is commonly 
done in the so-called Polyakov-loop-extended quark- 
meson models of three-color QCD [62-64]. To this end 
one introduces a constant temporal background gauge 
field Afj = Aod^o which is furthermore assumed to be in 
the Cartan subalgebra as in the Polyakov gauge, i.e., for 
the color SU(2) of QC2D simply given by Aq = T^2ao. 
This leads to the Polyakov loop variable 

0, = lTr,e'^^° = cos{J3ao), (15) 

to model a thermal expectation value of the color-traced 
Polyakov loop at an inverse temperature j8 = l/T", as an 
order parameter for the deconfinement transition at van- 
ishing chemical potential. With the covariant deriva- 
tive = djj - iSjjoAo replacing the ordinary one in 
the fermionic part of the QMD model in Eq. (13), this 
leads to a contribution of the form -iil/y^T^2aotl/ for 
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a Polyakov-loop-extended quark-meson-diquark model 
(PQMD) Lagrangian. 

It is often convenient to rewrite the fermionic part of 
the Lagrangian in yet another form, £g = 'J'g>S'~'*Pg, in 
terms of eight-component Nambu-Gorkov spinors *Pg = 
( t%'^ )' where tpr i^Pg) explicitly denote the red (green) 
color components of ij/ and if/'^ = CijF per flavor, and 

c-l _ / $-y''(iao+ii)+g(a-+iy^7i'f) gy^K \ ,^f-^ 

The PQMD model Lagrangian then becomes 

XPQMD = ^gS-'^g + 'Wpol(O) + h,d^O-f + h,d^7ff 

(17) 

where 'Z/poi(fl>) is the Polyakov-loop potential [42] 

which is commonly fitted to lattice results, but which 
can also be computed with functional methods [65, 66]. 
In contrast to the three-color case the Polyakov-loop po- 
tential is a function of one single real variable O here, 
even in the presence of a diquark condensate. 



2.2. Mesonic and baryonic fluctuations with the 
functional renormalization group 

The phases of QC2D and the competing dynamics of 
the fluctuations of the order parameters due to collective 
mesonic and baryonic excitations are most conveniently 
described within the framework of the functional renor- 
malization group (FRG) [67-73]. The central object in 
the Wilsonian RG flow equation as pioneered by Wet- 
terich [74], is the scale k dependent effective average 
action n [<I>], where <I> generically represents the set of 
all local fields in the theory. This scale-dependent ef- 
fective average action interpolates between the micro- 
scopic classical action at some ultraviolet (UV) cutoff 
scale k - K, aX which fluctuations of essentially all 
momentum modes are suppressed, and the effective ac- 
tion of the full quantum theory in the infrared (IR), for 
k ^ 0, which then includes all quantum and thermal 
fluctuations. The scale-dependence is described by the 
Wetterich flow equation. 



Tr{5,i?,(rf 



■RkT'] , (18) 



which involves a momentum- and scale-dependent reg- 
ulator R/i, whose precise form is not flxed but leaves a 
considerable flexibiUty. The role of the regulator is 



to suppress the fluctuations of modes with momenta be- 
low the renormalization scale k, and the flow equation 
is UV as well as IR finite. [O] are the second func- 
tional derivatives of the effective average action with re- 
spect to all the fields at scale k. The functional trace 
represents a one-loop integration typically evaluated in 
momentum space and includes the sum over all fields 
and their internal and space-time indices as well, with 
standard modifications for fermionic fields. It contains 
the full field and fe-dependent propagators of the regu- 
lated theory with cutoff R^, the inverse of [O] -1- R^. 
In order to solve the flow equation an initial microscopic 
action S = F^=a at some UV scale A has to be specified. 
Truncating the effective action to a specific form, the 
functional equation can be converted into a closed set of 
(integro-)differential equations, but will in general also 
introduce some regulator dependence in the flow. The 
choice of an optimized regulator minimizes this regu- 
lator dependence for physical observables. As bosonic 
(fermionic) regulators R^^b (Rkf) we choose 



RkAp) = (k^ - p'yeik^ - 



- 1 



e{k^-f). 



(19) 



which are three-momentum analogues of the opti- 
mized Litim regulators [75]. With this choice the 
three-momentum integration becomes trivial and the re- 
maining Matsubara sums can be evaluated analytically. 
Furthermore, this choice leaves the semilocal U(l)- 
symmetry of the Lagrangian unaffected, analogous to 
[50], where the chemical potential acts like the zero- 
component of an Abelian gauge field. 

At leading-order in a derivative expansion, all wave- 
function renormalization factors are neglected and only 
the scale-dependent effective potential Uk is taken into 
account. The ansatz for the effective average action then 
simply reads F^ = /(^^-^ XpqmdI^,,,^^^^.^^ ■ This means 
that in Xpqmd from Eq. (17) the meson/diquark poten- 
tial V0) of the 0(6) linear sigma model from Eq. (14) 
is replaced by Uk - ccr. The explicit symmetry breaking 
term -ca does not affect the flow and is thus not part of 
Uk but added after the RG evolution to the full effective 
potential again. At /i = 0, the scale -dependent Uk then 
only depends on the modulus of ^ = (cr, ji. Re A, Im A)^ . 
At non- vanishing chemical potential, however, we only 
have SO(4) X SO(2) symmetry and must therefore al- 
low it to depend on two invariants, i.e., Uk= Ukip^, cfi) 
where = cr^ + and (f = |Ap as before. For 

— > we recover the full extended SO(6) invariance, of 
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course, so that Uk then depends only on the combination 
(jp' — + (fi again. 

Working out the second functional derivatives of the 



effective action, one obtains in momentum space with 
constant fields, and coordinates such that cr - p, n - Q, 
Re A = f/, 1mA = 0, the inverse bosonic propagator 
r^B + ^M.with 



C = (20) 

'p^ + lUk,p ^ 

p^ + 2Uk,p 

p^+lUkj, 

' p' + 2Uk,p+^p'Uk„ 4pdUt,pd 

4pdUk^ p^ + 2Uk^+4d^Uk,dd-4fi^ -4m 

I Afipo p^ + 2Uk,d-4n^) 



where short-hand index notations for the derivatives of 
the potential with respect to the fields are defined as 



submatrix 



Uk,d = 



dUk 



,p = — , and Uk,pd = 



dp^ 



dpW 



etc. 



The inverse fermion-propagator can essentially be read 
off from Eq. (16). With the above replacements for 
constant fields it is diagonal in flavor and given by, 
Tf),+Rk,F, with 



n(2) 
k,F 



-'aoT+gp-TM 

-gr'd 



gy d 



v) ' 



(21) 



for each flavor With the regulators from Eqs. (19) the 
A:-dependent propagators are obtained upon inverting 
these expressions for Tfg^p+Rk^B/F- When these propa- 
gators are inserted into the Wetterich equation, Eq. (18), 
the integration over the spatial momentum components 
in the expUcit loop, with the three-dimensional opti- 
mized regulators in Eqs. (19), can be performed rather 
straightforwardly. For example, for the bosonic contri- 
bution to the flow, with bosonic Matsubara frequencies 
po - (jJn - 2nTn for periodic boundary conditions in 
the zero-component po of the loop momentum at finite 
temperature we have in the subspace of the three pions, 

(r<^(p) + RkAP)),. = {<^l + max}/, e] + 2Uk,p) 6ij 

with 1 , ... 3 for the upper-left 3x3 block of rfg 
in Eq. (20). In the subspace of the cr meson and the 
diquark-antidiquark pair, which mix when the diquark 
condensate d = |A| is non-zero, one analogously obtains 
a 3 X 3 matrix corresponding to the lower-right block 
in Eq. (20) with the same replacements, po a)„ and 
p^ — > max{/,A:^}. If we caU the inverse of this 3x3 



Aicon, max{f, k^]) = (r[^ + Rk,B)l\ , (22) 



the total bosonic contribution to the flow equation 
for the effective potential Uk (up to a 4-dimensional 
spacetime-volume factor V the same as that for the ef- 
fective action in this truncation), is then given by the 
following Matsubara sum 



dtUk,B = ^Tr { d,Rk,B [rf), + Rkb) 



(23) 



-2 Z^l 2n^ + 



-i-trA(a>„,A^)j 



6n^ 



oo _ 



The trace of A is the ratio of the sum of three 2x2 deter- 
minants of submatrices of (F^^^^^ + Rk,BiF)a-h o'^^r the 

full determinant of the 3x3 block (F^^^^^ + Rk,B/F)crA 
corresponding to the sigma and diquark-antidiquark di- 
rections in field space in Eq. (20). Therefore, the nu- 
merator of tr A is a quadratic polynomial in col and the 
denominator a cubic one. 



trA(c.„,^)= ^ (24) 



The coefficients of these polynomials are probably not 
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very illuminating. They are worked out to be [47], 

aro=3A:'' (25) 
+ 4fc2( - 4^/' + 2Uk4 + ^d^Uk4d + Uk,p + 2p^Uk,pp) 
+ 4(4/ + + 2UUd^Uk,dd + Uk,p + 2p''Uk,pp) 
- ^i^\Uk,d + d^Uk4d + Uk,p + 2p^Uk^pp) 
+ 2d\Uk,ddUk^ - 2p'ul^j + 2p'Uk.ddUk^p)) 
ai=6k^ + %Uk4 + ^(f'Uk4d + 4f/t^ + ^p^Uk^p 
/3o=(k'-4,/ + 2Uk4){k^ (26) 
+ 2k^{ - 2/|2 + Uk4 + 2(fUk4d + Uk,p + 2p^Uk,pp) 
+ 4( - 2i/Uk,p + Uk4Uk^ + 2(fUk4dUk^ 
- 4dVul^j + 2p\-2i/ + Uk4 + 2d^Uk4dWk,pp)) 
j8i = 3A:^ + 4k'{2Uk4 + 2d^Uk,dd + Uk,p + 2p'Uk,pp) 
+ 4(4/|4 + - An\Uk4 + d^Uk4d - Uk,p 
- 2p^Uk,pp) + 2Uk4{d^Uk4d + Uk^ + 2p^Uk,pp) 
+ 2d\Uk,ddUk,p - 2p^Ul^a + 2p^Uk4dUk,pp)) 
/32=3k^ + %t? + 4Uk4 + 4d^Uk4d + 2Uk,p + 4p^Uk,pp 

These coefficient functions a, and jS, depend on the 
renormaHzation scale, the chemical potential, the fields, 
and the derivatives of the potential. This mixing of the 
sigma with the diquark sector in the diquark conden- 
sation phase, where d - \A\ and baryon number is 
no-longer conserved, is what makes the equations some- 
what more complicated than usual. But the structure be- 
hind these lengthy expressions is actually not as bad as it 
might at first appear. With a partial-fraction decomposi- 
tion of tr A, for example, the Matsubara sum can still be 
computed analytically from the residue theorem in the 
standard way [76, 77]. With the roots of the polynomial 
in the denominator denoted as = -z?, i = 1,...,3, 
this yields for the bosonic flow 



k,B 



E" 

01 cothl-^ 
12;r2 {El \2T 



(27) 



ao 1 



coth 



(27) 



where E^ = yJW^Tw^p . 

The fermionic contribution to the flow can be worked 
out analogously. The derivative of the regulator Rkj in 
Eqs. (19) cuts off the spatial loop-momentum integra- 



tion at p - k , and for p < k the fermionic two-point 
function F^^^ + Rk^p is obtained from Eq. (21) upon re- 
placing -ip ■ f -i{p ■ ')'/|^) A: therein. Inverting the 



resulting fermion matrix, this then yields for the trace 
over Dirac (and flavor) indices for ^ < k^, 



ik^{Yfp+Rk,F 



\P\ 



(28) 



16r 



(pQ+aof-ix^+e+g\d^+p^) 
iipo + a^f + £f )((po + a^f + Ef) 



One factor of Nf = 2 hereby arises from the sum of 

the two flavors. The square of the denominator is the 

(2) 

determinant of the fermion matrix ^ -1- Rk,F per flavor, 
whose eigenvalues are given by the eight combinations 
of different signs of +j'(Po + ^^o) ± E^, with 



Et 



Jg^d^ + (ek ± /i)2 , and ek= Jk^ + g^p^ . 

(29) 

For [f- > k^ one would simply have to replace k^ by [f 
in these expressions again, as one would for F^^ with- 
out the regulator, e.g., for a mean-field analysis of the 
model, see Ref. [47]. These momenta are not needed 
for the optimized fermionic flow either, however, which 
then readily foUows to be given by, 

dtUk,F = -^Tvi^dfik,F{Tf^+Rk.Fy'} (30) 



(y„+ao)2-A/^+A:2+gV+p2) 
((y„ + ao)2 + Ef)(iv„ + oo)^ + Ef) 



where the fermionic Matsubara frequencies p[) = v„ = 
{2n+\)nT are used for antiperiodic boundary conditions 
in imaginary time. When their sum is evaluated, one 
finaUy obtains. 



dtUk,F = 



(31) 

{\-2Ng(E^/T,0)). 



Here, Nq(E IT ,Q>) are Polyakov-loop enhanced quark 
occupation numbers for QC2D, 

1 + (De^/^ 



NgiE/T,0). 



(32) 



1 -I- 20e'^'T + 

which reduce to the Fermi-Dirac and Bose-Einstein dis- 
tributions for O = 1 and O = -1, respectively. It is also 
common to use [42], 

sinhiE/T) 

<p^iE/T) = l-2Ar,(£/r,(D) = '' . (33) 

0) + cosh(£/r) 

which satisfies 

' tanh(x/2) = 1 - 2p/(x) , 0) = 1 
ip<t>{x) - I tanh(x) = 1 - 2pj(2x) , d) = 
coth(x/2) = 1 + 2pfc(x) , 3) = -1 

(34) 
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where = — — - , and pi,(x) = — — - . 

The SU(2) center sector with O = -1 corresponds to 
using periodic quarks in QC2D and leads to a bosonic 
excitation spectrum. It is expUcit in Eq. (30) that noth- 
ing changes, if one changes to periodic quarks and 
changes the background according to such a center flip 
O — > -O at the same time. The first amounts to replac- 
ing v„ oj„ and the second to uq ^ uq + nT which 
together leave the fermionic flow invariant.^ This is a 
manifestation of the so-called Roberge-Weiss symme- 
try [78]. For C> = 0, in a center- symmetric background 
with flo - nT 12 mod n, the thermal excitation energies 
are twice the quark energies which is occasionally inter- 
preted as modeling confinement. 

The full and final flow equation for the effective po- 
tential in the PQMD model for QC2D [47] in flie lead- 
ing order derivative expansion is given by the sum of the 



bosonic and the fermionic flows in Eqs. (27) and (31), 

5,f/* = d,lJk,B + d,Uk,F . (35) 

A significant complexity in this flow equation for the ef- 
fective potential U (p, d) is introduced by the presence of 
two fields p and d corresponding to the standard chiral 
quark condensate and the diquark condensate, respec- 
tively, which have a mutual influence on one another 
with physically important impUcations. 

2.5. 50(6) symmetric flow 

In the normal hadronic phase without diquark con- 
densation, we may set J = |A| = to obtain an ex- 
plicitly SO(6) symmetric flow for f/j(0). If we set 
IJ]^^^ = Uk^ = Uk4, Eq. (35) then reduces to a more 
familiar looking form. 



{ 3 



£7 



IT 



1 



\ — coth + — coth + ™ coth 



2T 



1 



\ 2T 



- ^ ( 1 - Af, (e, - //; r, O) - N^iek + fi; T, O) )| , 



(36) 



with equal single-particle energies for mesons and di- 
quarks, 

' (37) 



while for the sigma meson one has 

£^ = + 2Uk4 + 4<*2f/*.# . (38) 



and for the quarks - yjk^ + g^cf^ as above. Except 
for the change in the number of active degrees of free- 
dom contributing to this flow, and the isospin-like chem- 
ical potential coupling to one pseudo-Goldstone boson 
pair, the SO(6) synametric flow equation here is entirely 

^In Eq. (30) we can actually continuously change v„ -» v„ + c and 
ao -> ao - c without effect. It might appear that we can thus change 
the quarks' boundary conditions by a general U(l) phase as we could 
in QED with unbroken displacement symmetry. This is not the case, 
however, because v„ — > v„ +c; would rotate the boundary conditions of 
the red quarks in the same way as those of the green antiquarks, i.e., 
those of the red ones in the opposite direction of those of the green 
ones. Thus this U(l) rotates about the T^-direction of SU(2) in color 
and there is no contradiction. 



analogous to the one of the PQM model for QCD with 
three colors, see e.g., [64, 79, 80]. In the PQM model 
for QCD with isospin chemical potential, however, one 
must allow for pion condensation and then arrives at a 
flow equation for two competing fields [57] analogous 
to Eq. (35), with Eqs. (27) and (31). 

A further simplification occurs for // = and in a 
trivial gauge-field background, ao = corresponding to 
0=1, 



k^ j 5 
\27fi \El 



coth 



ZL 
2T 



1 



coth 



2T , 



(39) 



Because // = the diquarks are now fully degenerate 
with the pions which leaves us with the Nc = 2 ana- 
logue of the familiar three-color QM model flow equa- 
tion [81, 82] except that there are now five pseudo- 
Goldstone bosons instead of the usual three pions. 

As in the studies of 0(4) universality and scaling in 
the three-color QM model [9-11, 83, 84], one can use 
this flow equation to analogously check the synametry 
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Figure 4: The chiral condensate in the QMD model for QCtD over 
temperature T from the SO(6) symmetric RG flow for various quark 
masses con'esponding to pion masses = 138 MeV, 84 MeV and 
in the chiral limit, relative to the 7 = condensate at = 138 MeV. 

breaking pattern discussed in Section 2.1, and the corre- 
sponding magnetic scaling, by computing the two criti- 
cal exponents p and 6 as mentioned in the Introduction. 
For/v = m - 0, the SU(4) ^ SO(6) dynamically breaks 
down to Sp(2) ^ SO(5) so that we expect a finite tem- 
perature phase transition in the three-dimensional 0(6) 
universality class. The critical exponent /3 can be ex- 
tracted from the dependence of the chiral condensate on 
the reduced temperature t - (T - Tc)/Tc in the chiral 
limit, whereas the exponent S governs the dependence 
of the chiral condensate at on the quark mass irig or 
correspondingly on the explicit symmetry-breaking pa- 
rameter c, 

{qqh ~ i-tf, {qq)T, ~ c'l' . (40) 

With the usual two-exponent scaling all other critical 
exponents are then obtained from these two. A careful 
analysis of the behavior of the quark condensate from 
solutions to the \d flow equation (39) such as those 
shown in Fig. 4 yields critical exponents [3 - 0.4318(4) 
and 6 - 5.08(8) [47]. Literature values for these expo- 
nents obtained from Monte-Carlo simulations are given 
by = 0.425(2) and 5 = 4.77(2) [85]. At leading or- 
der in the derivative expansion one would not expect 
to reproduce these values exactly, however The more 
appropriate benchmark here should be the functional 
renormalization group result for the 0(6) model in the 
local potential approximation [86]. In absence of wave- 
function renormalizations there is no anomalous dimen- 
sion for the fields and their critical exponent therefore 
vanishes, 77 = 0. Then the hyperscaling relations, 

P-'.id-2.,), (41) 
a - 2 + J] 2 
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Figure 5: The degenerate pion and diquark masses in comparison with 
that of the sigma meson over temperature T from the S0(6) symmetric 
RG flow for the same model parameters as in Fig. 4. 

immediately entail that 6 = 5 and yS = v/2 in li = 3 
dimensions. With the correlation-length critical expo- 
nent V = 0.863076 from Ref. [86] this corresponds to 
/3 = 0.4315, and both these values are in good agree- 
ment {i.e., within the errors) with the two from Eq. (39), 
yS = 0.4318(4) and 6 = 5.08(8) as quoted above [47]. 

The corresponding temperature dependent screening 
masses as extracted from the resulting effective poten- 
tial also look qualitatively exactly the same as in the 
QM model and are shown in Fig. 5 for the same set of 
zero-temperature pion masses as in Fig. 4. 

If we dismiss for the moment the effects of a finite 
density of the bosonic baryons of QC2D, and the cor- 
responding fluctuations due to the collective baryonic 
excitations at finite density, we may go ahead and solve 
the one-dimensional SO(6) symmetric flow equation in 
Eq. (36) at finite quark chemical potential // and finite T. 
The resulting quark condensate is shown in the three- 
dimensional plot in Fig. 6. The corresponding phase 
diagram is shown in Fig. 7. 

The most important aspect of these results is that they 
are again qualitatively exactly as in the quark-meson 
model for QCD. There is a low-temperature P' -order 
transition line with an endpoint at around // « 2.5 m„, 
and a chiral crossover which is here indicated by the 
dashed lines marking the half-value of the vacuum con- 
densate. As also known from studies of the quark- 
meson model [82], the inclusion of fluctuations in the 
order parameter to account for collective excitations 
with the ERG leads to the capacity to describe a low 
temperature phase transition to bound quark matter. As 
shown in Figs. 6 and 7, chiral symmetry is not fully re- 
stored but only jumps at the transition and gradually de- 
creases further into the quark-matter phase. This is con- 
ceptually different from the restoration of chiral sym- 
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Figure 6: The chiral condensate from the SO(6) symmetric RG flow 
over quark chemical potential and temperature. The /j = plane cor- 
responds to the data shown in Fig. 4 for the = 138 MeV pion mass. 



metry at the first-order line in mean-field studies which 
would lead at the quark level to the analogue of Lee- 
Wick matter as in the chiral Walecka models [87], but 
which is a mean-field artifact. 

Another nice feature of the SO(6) symmetric FRG re- 
sult at finite fi is that the zero-temperature quark con- 
densate remains precisely constant until jj reaches the 
critical value for the transition to bound quark matter as 
seen in Fig. 6. There is no Silver-Blaze problem in the 
SO(6) symmetric flow from Eq. (36). The vacuum does 
not change at all until the quantum phase transition to 
quark matter is reached. 

These results are all nice and well to illustrate general 
features of quark-meson models, but they fail of course 
to describe an essential part of the dynamics which are 
the effects of finite baryon density and collective bary- 
onic excitations. In QC2D this means to correctly de- 
scribe the diquark condensation phase with superfluidity 
of the bosonic baryons and the BEC-BCS crossover in 
the cold 'quarkionic' two-color quantum gas."^ The ad- 
vantage here as compared to the three-color case is that 
we can understand these effects quite well and study in 
detail what we need to do to include them. Moreover, 
it is not only much more straightforward to include the 
baryonic degrees of freedom in form of the diquarks of 
QC2D than the true fermionic baryons in the real world 
with three colors, but it can also be considered a valu- 
able warm-up exercise for a successful quark-diquark 
description of baryonic degrees of freedom in QCD. 



QMDMF ■ 
QMDRG ■ 



1.5 



"^The emphasis in 'quarkionic' is on replacing the fermions in the 
ultracold fermionic gases by the colored quarks which are unphysical. 
If it is on baryonic excitations as the physical ones in confined quark 
matter, quarkyonic may be also appropriate. 



Figure 7: QMD model phase diagram without diquark condensation 
and baryonic fluctuations. Shown are results from [47] for the S0(6) 
symmetric flow (RG) and a typical mean-field diagram including vac- 
uum contributions (MF). 



2.4. Diquark condensation 

At zero temperature, the vacuum and the mass 
spectrum must remain independent of the (relativistic) 
chemical potential ju until the quantum phase transition 
is reached. To demonstrate that this phenomenologi- 
cal feature holds in the functional integral representation 
of the grand potential has been called the Silver Blaze 
problem after the story by Arthur Conan Doyle in which 
the dog did not bark when the horse named Silver Blaze 
was abducted [52]. This "curious incident of the dog in 
the night-time" doing nothing led Sherlock Holmes to 
solve the case. The solution to the Silver Blaze problem 
for QCD with isospin chemical potential [52] equally 
holds for QC2D with baryon chemical potential. How- 
ever, the relevant quantum phase transition in QC2D is 
not the QM model transition to quark matter, but that of 
diquark condensation at /i = m„/2 with the spontaneous 
breaking of baryon number and baryon superfluidity. 

Diquark condensation at jj = m„/2 thus occurs way 
before any quark matter or chiral transition at around 
yu ^ 2.5 m„ as seen in the essentially purely mesonic 
model above. This is well known from chiral effec- 
tive field theory [21], lattice simulations [31] and mean- 
field studies of the NJL [40] or die PNJL model [42] as 
mentioned at the beginning of this section. In order to 
describe it within the FRG one needs to solve the full 
flow equation (35) with including the competing fluctu- 
ations in the two directions of field space corresponding 
to the chiral and diquark condensates, respectively, both 
at the same time. This was done numerically on a two- 
dimensional grid in field space in Ref. [47]. 

Of course, it can also be described qualitatively al- 
ready at mean-field level in the QMD model [47]. The 
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resulting meson and diquark mass spectrum is shown 
in Fig. 8. Just as the dichotomy between chiral sym- 
metry breaking and the Goldstone theorem is manifest 
within the mean-field/RPA framework or in the rain- 
bow/ladder truncation to Dyson-Schwinger equations, 
the Silver Blaze problem is absent here as well. The 
meson masses at zero temperature remain constant at fi- 
nite jj below the transition. It can be shown analytically 
that the diquark and antidiquark masses are given ex- 
actly by - m„ + 2jj there, and the independence of 
the fermion determinant on fi for yU < m„/2 follows. 

It is important to stress that at the mean-field level, 
this exact feature of the theory holds only for the pole 
masses in the meson and diquark propagators as ob- 
tained from the random-phase approximation. This is 
true in the (P)QMD model as it is in the (P)NJL model. 
The differences between the two are not significant and 
only quantitative in nature in this approximation. The 
corresponding pole masses are defined as the zeroes of 
the bosonic two-point function, 

detr^^'(;?) = 0, for - ^ mj , k^l,...6. (42) 

Defining 

rf(p)^r%(p) + TlB(p), (43) 

(2) 

where F^^j and 11^ are the tree-level and vacuum- 
polarization contributions, respectively, and with the 
bosonic potential V{(f) from the linear-sigma model, 
Eq. (14), in the normal phase without diquark condensa- 
tion the pole-mass conditions (42) are then simply given 
by the solutions of the following equations for a>, 

: Li/ = -m^ + 3Acr^ + Ylg-i'^, T) , 

m„ : (!?■ - -m^ + Acr^ + I\-„{(x>, T) , 

mA± : (w ± l^if = -m^ + Acr^ + nA±(w, T) . 

(44) 

(21 

Up to the chemical potential entries in F^ Jj the bosonic 
two-point function is diagonal in the normal phase with 
d - Q, and we have introduced n,(tij, T) for the six di- 
agonal entries of YIb(p) with p - (-ico, (J) at finite T. 

As with chiral symmetry breaking and the occurrence 
of massless pions in the chiral limit, one then readily 
verifies that the gap equation for the diquark condensate 
along the bifurcation line p.c{T), which defines the on- 
set of U(1)b breaking and diquark condensation in the 
(r, //)-plane, leads to the condition, 

( - + Acr^ - 4// + Ha* (0, T))\ = , (45) 

because the corresponding derivative of the effective po- 
tential is equal to nA+(0, T) = nA-(0, T). This im- 




IX [m„] 

Figure 8: Pole-mass spectrum at T = 0: mean-field/RPA QMD model 
results from Ref. [47]. For comparison also shown: without the efi'ect 
of diquark/sigma-meson mixing in the superfluid phase, linear sigma 
model and^PT results. 

plies that a solution with diquark-mass zero exists in 
Eqs. (44) along the same line pdT) in the (T,p)-plane.^ 
One can further verify analytically that Ha* (oj, T) - 
Tl„(cj, T) for yt/ = as it must from SO(5) symmetry, 
and, more significantly, that at zero temperature, 

Ha* (w, 0) = n^(co ±2/^,0) . (46) 

Therefore, the T - solutions to Eqs. (44) obey 

ota± — m„ ± 2ii . (47) 

Together with the y[/-independence of the pion and 
sigma masses, this guarantees that there is no Silver 
Blaze problem at mean-field/RPA level as seen in Fig. 8, 
and that the zero-temperature quantum phase transition 
occurs atp - fidO) - m„/2. 

In contrast, in quark-meson models one usually as- 
signs the physical masses of the mesons to the so-called 
screening masses to fix the model parameters. These are 
obtained as the eigenvalues of the Hessian of the effec- 
tive potential in the leading-order derivative expansion. 
They agree with the eigenvalues of 

4'^(0) = F(,'ji(0) + nB(0), (48) 

and with the correct analyticity properties [76, 77] in the 
continuation from discrete Matsubara to continuous real 
frequencies this equivalence holds at all temperatures. 



^^This Goldstone mode persists everywhere into the diquark con- 
densation phase but the simple pole-mass conditions (44) are no- 
longer valid in this form there. 
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Figure 9: Zero temperature condensates from full flow compared to 
mean-field and^PT results, and lattice data (from [31]). 

However, it is straightforward to verify [47] that the 
screening masses have an enormous Silver Blaze prob- 
lem. Their behavior with chemical potential is utterly 
unphysical. That of the sigma meson does not remain 
independent of yt^ at T = in the normal phase as it 
should. The diquark and antidiquark screening masses 
remain degenerate throughout the normal phase at all 
T, and they both vanish at fidT). However, ndT) is 
not half the screening mass of the pion but half its pole 
mass. One would have thought that the difference is 
small, as we only have to extrapolate from = to 
= -ml in the momentum dependence of the two- 
point function to get from one to the other, but this is 
not the case. In the QMD model of QC2D, where the 
the quantum phase transition at fidO) - m„l2 defines the 
physical pion mass, one observes that the corresponding 
screening mass is about 30% too large [47]. And this is 
true both in mean-field/RPA and with the effective po- 
tential obtained from the full flow equation (35) with 
mesonic and baryonic fluctuations. One thus needs to 
compute the pion pole-mass in those models to fix the 
parameters to more realistic values. At the mean-field 
level an RPA computation will do. A pole mass calcu- 
lation from a flow equation for the fully momentum de- 
pendent r^^(p) in a truncation consistent with that for 
the effective potential indeed shows that this is much 
closer to m„ = l/j^ than the screening mass from the 
effective potential in the FRG computation as well [47]. 

In the diquark-condensation phase the sigma meson 
mixes with the two diquark modes, i.e., the respective 
masses have to be determined from the zeroes of the de- 
terminant of the corresponding 3x3 submatrix in . 
As in the NJL model [42, 55] one can verify further ex- 
act results from the mass formulas at T = 0. Also in 
the QMD model at mean-field/RPA level the in-medium 



Figure 10: The chiral and diquark condensates from the full 2d RG 
flow over quark chemical potential and temperature. The p = Q plane 
corresponds to the data shown in Fig. 4 for the m,r = 138 MeV pion 
mass, the 2d grid solution with explicit diquark fluctuations agrees 
with that of the SO(6) symmetric flow as shown in Fig. 6 there. 

pion pole-mass is equal to m„ - 2fi above the onset 
of diquark condensation at 2/z = m,t,o- Moreover, one 
verifies explicitly that one of the three modes in the di- 
quark/sigma sector remains exactly massless in the su- 
perfluid phase, also at finite temperature. This is of 
course the Goldstone boson corresponding to the spon- 
taneously broken U(1)b baryon number. Another one 
becomes degenerate with the pions for large values of 
the chemical potential, eventually, reflecting the restora- 
tion of chiral symmetry. They combine into an SO(4) 
multiplet as the chiral condensate vanishes for large fi. 

Fig. 8 also shows the corresponding mass formulas 
from the leading order chiral Lagrangian and the linear 
sigma model along with results of an RPA mass cal- 
culation where the mixing terms in the sigma/diquark 
sector were neglected. In those results all off-diagonal 

(2) 

terms in F^ other than the chemical potential entries 

(2) 

in F^ jj were set to zero by hand to demonstrate that it 
is these terms which lead to an avoided crossing in the 
sigma/diquark channel. In particular, without this mix- 
ing, the masses show more clearly that it really is the 
original sigma which becomes degenerate with the pion 
at high density as required by chiral symmetry. The 
transition marks the BEC-BCS crossover in which the 
massive diquark mode changes from its original would- 
be-Goldstone nature (shifted by 2/i) to the heavy Higgs 
mode in the BCS limit at large jj.. 

From the discussion above, it is clear that it is best 
to fix the pion mass to the zero temperature quantum 
phase transition in the full FRG calculation of the ef- 
fective potential from Eqs. (35), (27) and (31). A cor- 
responding pole mass can be calculated from the flow 
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Figure 1 1 : Phase diagram from RG flow witli collective baryonic fluc- 
tuations (and no chiral 1"^' order transition/critical endpoint) compared 
to the SO(6) symmetric flow in Fig. 7 without diquark condensation. 



Figure 12: Comparison of QMD phase diagrams from MF and RG 
calculations, including lines with gcr = p in the superfluid phase to 
indicate the BEC-BCS crossover. 



of the bosonic propagator but it is not needed once we 
know they must agree. The resulting zero temperature 
chiral and diquark condensates are shown in Fig. 9. The 
linear sigma model expressions for the T = conden- 
sates cr = {qq} and A = {qq} are [43] 

O" _ J 1 for fJ. < iJ^c 

^"1 ^\-^,+2fr\ for ti>ij, ' 

where x - 2fi/m„ and y = ma-/m„. The only difference 
between these and the leading order ;^'PT result [21], i.e. 
in the non-linear sigma model, is the y-dependent term 
in the diquark condensate which reduces to the;t'PT for- 
mula for y — > CO. In this limit the vacuum realignment 
from (^g')-like to and (<7<7)-like is described by a simple 
rotation with constant {qq}^ + {qq}^. 

2.5. Phase diagram of the QMD model for QC2D 

The solution to the full flow Eq. (35), with Eqs. (27) 
and (31) on a two-dimensional grid in field space can 
be extended into the entire (r,jt/)-plane to demonstrate 
the effect of baryon density on the chiral condensate 
in QC2D [47]. For vanishing gauge-field background, 
O = 1, the result is summarized in Fig. 10. This figure, 
from the full two-dimensional flow for an eff'ective po- 
tential with the reduced SO(4)xSO(2) symmetry, should 
be compared to Fig. 6 from the flow for the SO(6) sym- 
metric effective potential, to appreciate the influence of 
finite baryon density. 



For small chemical potential in the normal phase the 
results agree very well. Considering that the full so- 
lution does include independent mesonic and baryonic 
fluctuations also in this regime, albeit the latter around 
A = 0, this is perhaps not quite so obvious as one might 
think at first. It confirms, however, that a purely mesonic 
model can produce reliable results so long as the baryon 
density remains zero or sufficiently small, i.e., the part 
of the phase diagram relevant for the mesonic freeze- 
out. Those results are unambiguously determined by 
the chiral symmetry breaking pattern, here in the 0(6) 
universality class. Allowing additional interactions with 
lower symmetry has no effect on the flow in this regime. 

Once the quark-chemical potential approaches half 
the baryon mass, corresponding to ms/Nc, however, the 
rapidly increasing baryon density equally rapidly sup- 
presses the chiral condensate. With the proper inclu- 
sion of the collective baryonic excitations, there is no 
trace left of the chiral first-order transition and the crit- 
ical endpoint of the purely mesonic model. At least at 
zero temperature, the baryon density is an order param- 
eter for N f - Nc - 2 as well, and it rapidly increases at 
the finite temperature transition line separating the nor- 
mal from the superfluid phase which should give rise to 
the two-color analogue of the baryonic freeze-out. 

The corresponding phase diagrams from the two RG 
solutions with and without finite baryon density are 
compared in Fig. 11. The dashed Unes again denote the 
half-value of the chiral condensate in the vacuum as an 
indication of the chiral crossover line. 

The onset of diquark condensation and superfluid- 
ity of these bosonic baryons at low temperatures also 
marks the line at which the the residual SU(2)i/ flavor 
symmetry starts changing in nature from the approxi- 
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mate Sp(2) symmetry in the normal phase to becoming 
the approximate SU(2)l x SU(2)r quasi-restored chiral 
symmetry. Because they are both exphcitly broken and 
only approximate symmetries, this vacuum realignment 
is a crossover. The quark mass with large chiral conden- 
sate in the normal phase starts out as the predominantly 
spontaneously generated Dirac mass, and the bosonic 
baryons undergo Bose-Einstein condensation in form of 
a dilute gas of strongly bound diquarks at the onset of 
diquark superfluidity. As their density increases, the un- 
derlying quark mass rotates into a spontaneous Majo- 
rana mass leading to a BCS-like pairing. This is the rel- 
ativistic analogue in two-color QCD of the BEC-BCS 
crossover in ultracold fermionic quantum gases. It is 
indicated in Fig. 12 by additional dashed lines in the 
superfluid phase where the quarks' Dirac-mass equals 
their chemical potential. For a more comprehensive dis- 
cussion within the NJL model, see Ref. [46]. 

In Figure 12 the phase diagram from the full RG so- 
lution is compared to a QMD model mean-field result 
from Ref. [47]. As temperature increases, the line of the 
diquark-condensation phase transition, which should be 
in the 0(2) universality class, in the FRG solution with 
fluctuations differs more and more from the mean-field 
result. The first-order transition line at high temper- 
atures in the mean-field approximation is washed out 
by the fluctuations, and the associated tricritical point 
which was as also predicted from next-to-leading order 
XPT [23] turns out to be a mean-field artifact. 

To summarize, the particular advantages of using two 
instead of the usual three colors here were twofold: 

- to prepare for the inclusion of diquark correlations 
and baryonic degrees of freedom in a covariant quark- 
diquark description, e.g., by a corresponding quark- 
meson-baryon model for real QCD as a next step. 

- to be able to test non-perturbative functional methods 
and models against exact results and lattice simulations. 

The main results from the tests discussed in this sec- 
tion include: the verification of 0(6) scaling at yu = 0; 
the demonstration of the relevance of pole masses to 
correctly describe the zero temperature quantum phase 
transition of two-color QCD in the FRG framework, and 
the failure of the usual screening masses to be capable of 
that; and finally but most importantly, the non-existence 
of a chiral first-order transition and critical endpoint at 
finite baryon density. 

3. Universal aspects of deconflnement 

The phases of gauge theories can be classified accord- 
ing to the behavior of the flux of a static fundamental 
charge. In a Coulomb phase, the electric flux spreads 



. Boundary conditions: j^,, .j Hooft s 

anti {C-)periodc twisted b.c.'s twisted C-periodic b.c.'s 




Coulomb Confinement Higgs 

V y 

mass gap 

Figure 13: The phases of gauge theories classified by the flux of a 
single static test charge alone in the world. 

out and extends to infinity because of Gauss' law. The 
Coulomb cloud of the charge is isotropic and of long- 
range nature because the photon is massless. With a 
mass gap, there are no such long-range Coulomb fields. 
Then the distinction between Higgs and confinement 
phases depends on whether charge is conserved or not. 
If it is, there is a kind of Gauss law and the flux of the 
test charge must go somewhere. Because of the mass 
gap, it will get squeezed into a string to minimize the 
cost of free energy. In a Higgs phase, on the other hand, 
charge is not conserved, there is no Gauss law, the con- 
densate screens the charge and its flux lines peter out. 
This rough classification is sketched in Fig. 13. In the 
language of local quantum field theory, in terms of local 
field systems to measure the colored quantum numbers 
of the elementary degrees of freedom but not necessar- 
ily with physical asymptotic states associated to these 
fields (which would be called field-particle duality), this 
is described by the Kugo-Ojima criterion. This criterion 
specifies the same two conditions for confinement in the 
local field theory language: (i) there must be a mass gap 
to avoid the analogue of the Coulomb cloud in QED, but 
(ii) all global gauge charges must be conserved to have 
a Gauss law and to avoid the Higgs mechanism. 

When the charges of global gauge invariance are un- 
broken one derives more generally that every gauge- 
invariant localized state is a singlet under these un- 
broken global gauge charges. Thus, without (elec- 
tric) Higgs mechanism, QED and QCD have in com- 
mon that any localized physical state must be charge- 
less/colorless. The extension to all physical states is 
possible only with a mass gap, however Without that, 
as in QED, non-local charged states which are gauge- 
invariant can arise as limits of local ones which are not. 
The Hilbert space decomposes into so-called superse- 
lection sectors of physical states with different charges. 
A sector of (total) charge one differs from the neutral 
sector by the presence of one charged particle, and the 
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energy difference between the ground states of both sec- 
tors defines its mass. With a mass gap in QCD, on the 
other hand, color-electric charge superselection sectors 
cannot arise: every gauge-invariant state can be approxi- 
mated by gauge-invariant localized ones (which are col- 
orless). One concludes that every gauge-invariant state 
must also be a color singlet. 

To prepare a charged sector, one starts in a finite 
spatial volume with suitable boundary conditions. The 
free energy difference between the charged and the neu- 
tral sector then approaches a constant finite value, zero, 
or it tends to infinity in the infinite volume limit, in a 
Coulomb, a Higgs, or a confined phase, respectively. 

This is all well established and tested in pure gauge 
theories with infinitely heavy test charges but without 
dynamical (color-)charged particles. Beyond that, in 
particular at finite temperature it is not so clear, how- 
ever. A mass gap can change gradually so the high tem- 
perature deconfined and the low temperature confined 
phases can be analytically connected as in full QCD 
with a crossover. But even at zero temperature where 
there should be a clear distinction between differently 
charged superselection sectors, there is not. The meth- 
ods to fix the total color charge which will be described 
below only work for the pure SU(A') gauge theory with- 
out dynamical quarks, at least when QCD is studied as 
an isolated theory. An alternative where these methods 
could be applied might be to go beyond that and, in par- 
ticular, include the quarks' fractional electric charges 
as in the Standard Model. This possibility will be dis- 
cussed briefly in Subsection 3.5. 

In QCD alone, a non-zero total fundamental color 
charge cannot be screened completely without Higgs 
mechanism when dynamical quarks are present either, 
but with string breaking it should become a bound- 
ary effect and the energy difference to the neutral sec- 
tor should reduce to the corresponding string-breaking 
scale in the thermodynamic hmit. Moreover, one might 
expect that the Kugo-Ojima criterion should predict 
symmetry breaking and thus a phase boundary between 
confined and Higgs phases, for example, but there are 
counter examples on the lattice for that as well [88]. 

In pure gauge theories where the methods to fix 
differently charged sectors work well, we often also 
have a duaUty between magnetic and electric sectors. 
For example, antiperiodic (spatial) boundary conditions 
in the compact U(l) lattice gauge theory, where they 
are equivalent to C-periodic ones, can be used to pre- 
pare magnetic sectors and compute the corresponding 
monopole mass [89] to distinguish its Coulomb and 
confined phases. The difference in free energy of the an- 
tiperiodic versus the periodic ensemble thereby tends to 



zero or a finite value for the confining (magnetic Higgs) 
or Coulomb phases, respectively. By duality, via the 
Z gauge theory with integer valued link variables, the 
(electric) Higgs vs. Coulomb phases of the non-compact 
Abelian Higgs model follow an analogous pattern. 

In SU(A^) gauge theories the analogue of antiperi- 
odic boundary conditions, which would not in general 
be gauge invariant there, are 't Hooft's twisted bound- 
ary conditions to fix the total Zff center-flux through the 
planes of the spacetime torus by fixing the total num- 
ber of center vortices mod through these planes [90]. 
Dual to those are the corresponding electric flux sec- 
tors which describe the free energy Fg{T, L) of the flux 
of a static fundamental charge in a l/T x box with 
boundary conditions to mimic the presence of a mirror 
charge in the neighboring volume along the direction of 
the flux. This is a bit more complicated than the sim- 
ple anti/C-periodic boundary conditions in the Abelian 
case but it has the analogous interpretation in terms of 
mirror charges. It is reviewed in Subsec. 3.2. Moreover, 
combinations of C-periodic and twisted boundary con- 
ditions can be used to calculate the mass of the 't Hooft- 
Polyakov monopole in SU(A^) with adjoint Higgs at 
least for even [91]. 

3. 1. Center vortices and spin interfaces 

The finite temperature deconfinement transition in 
SU(A^) gauge theories \n d + I dimensions is very well 
understood in terms of the spontaneous breakdown of 
their global Z^r center symmetry [92]. This symmetry 
is faithfully represented by the fundamental Polyakov 
loops P{x) which five on the d spatial dimensions and 
describe static fundamental charges. Under the global 
Zjv center symmetry they transform like spins Sj in a c/ 
dimensional ^-state Potts model with q = N and Hamil- 
tonian [93], 

9{^-jYj^s,sj-hY,^s>,o, Si = 0,l,...q-1 , 

(i.j) ' 

(50) 

with nearest neighbor couphng /. A non-zero external 
field H, inversely related to the quark mass niq, may be 
included to mimic the leading effect of heavy dynamical 
quarks. When l/rriq = 0, the Polyakov loop develops a 
non-zero expectation value only in the deconfined, Z^y- 
broken phase, while the expectation value of van- 
ishes in the disordered, confined phase much Uke the 
spontaneous magnetization in the spin model. 

This is well described in terms of spacelike center 
vortices which play the role of spin interfaces. They 
separate regions where the Polyakov loop differs by 
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{d-H) dimensional SU{N,) 




J spatial dimensions d 
self-dual for all q 



Figure 14: Separation line between P' and 2 order phase transitions 
in the ^-state Potts models, from [93], and pure SU(A') gauge theories. 

a phase z e Z^, and their proliferation disorders the 
Polyakov loop and leads to confinement. The suppres- 
sion of spatial center vortices at high temperatures coin- 
cides with the ordering of the Polyakov loop, and their 
free energy offers an elegant order parameter for the 
transition [94]. 

In order to be able to apply the powerful tools of uni- 
versality and scaling near a critical point, here we are 
particularly interested in cases where this transition is 
of 2"'' order For pure QCD with N = 3 colors in 3 + 1 
dimensions this is of course not the case. The transi- 
tion is first order, but only just. If we either reduce the 
number of colors to = 2 or the dimensions to 2 4- 1, 
or both, we obtain Yang-Mills theories with a second 
order deconfinement transition within the universality 
class of ^-state Potts models with q = 2 (Ising) in li = 3 
dimensions or with q - 2 and 3 in li = 2 [95]. The 
q - 4 Potts model in 2 dimensions is interesting be- 
cause it is known [96] to have a 2"'' order transition and 
to fall precisely on the separatrix qdd) with respect to 
the order of the transition in the (q, ii)-plane as shown 
in Fig. 14, where only P' order transitions occur for 
q > qdd), i.e., qc{2) - 4. The corresponding SU(4) 
gauge theory in 2-1-1 dimensions has been studied for 
example in [97-99]. The conclusion in [98] was that 
the transition is weakly I'*' order, unlike the q - A Potts 
case. It nevertheless seems that there is a rather wide 
range of intermediate length scales where at least ap- 
proximate Potts scaling can be observed [97, 100]. One 
might then like to understand, for example, why among 
the wider class of Z4-symmetric Ashkin-Teller models 
with continuously varying critical exponents, it is the 
standard q - A Potts scaling that is relevant here, and 
whether this can be derived from an effective Polyakov- 
loop model analogous to the known cases with N - 2 
and 3, along the fines of Refs. [101-103]. 

One general property of the 2d Potts models is that 



they are selfdual for all q. We will see below, how this 
selfduality is reflected in a duality between the spacelike 
center vortices and the confining electric fluxes of the 
gauge theory. 

In the pure SU(A^) gauge theories 't Hooft's twisted 
boundary conditions fix the total Z/vi-valued center flux, 
i.e., the total number modulo of center vortices in 
the various planes of the Euclidean spacetime box [90]. 
At finite temperature T, the free energy differences of 
these twisted ensembles and the periodic one define 
the corresponding center- vortex free energies. If they 
tend to zero in the thermodynamic limit, then these vor- 
tices condense which leads to an area law for those 
Wilson loops that feel their disordering phases in the 
corresponding plane of the twist. In a l/T x L'' box, 
we have to distinguish between temporal and magnetic 
twist. The latter is defined in purely spatial planes and 
corresponds to the Z^f-valued magnetic flux m of static 
center monopoles, with the direction of m perpendicular 
to the spatial plane of the twist. Their free energy tends 
to zero with L — > oo as tx^{-crs{T)L^] at all T [104], 
corresponding to the area law for spatial Wilson loops 
with spatial string tension crdj). 

Relevant for the deconfinement transition are only 
the temporal twists. These are characterized by d- 
dimensional vectors k of integers mod A^, i.e., with 
components ki e Z^i representing the twist in the tem- 
poral plane of orientation (0, /), with total center flux 
exp(2m ki/N) through that plane. See Fig. 15 for an 
illustration of such a vortex ensemble in SU(2). We 
denote the partition functions of the various ensembles 
with temporally twisted boundary conditions by Zk{k). 
The corresponding center-vortex free energies per T are 
then given by 

F^(jl) = _ inRdk) = - In (Z,(^)/Z,(0)) , (51) 
ki^O,l,...N -I , 

where Zi(0) stands for the periodic ensemble. 

Intuitively, center vortices can lower their free energy 
by spreading out. As temperature is increased, how- 
ever, the temporal ones can no-longer spread arbitrarily 
but get squeezed more and more until the phase tran- 
sition is reached above which they are completely sup- 
pressed. In the vicinity of a second order deconfinement 
transition these vortex free energies show the universal 
behavior of interfaces in the respective tZ-dimensional 
Potts model. Interfaces in the spin models are typically 
introduced as frustrations along which the coupling of 
adjacent spins favors cyclically shifted spin states rather 
than parallel ones for the usual ferromagnetic couplings 
y > 0. They form d - I dimensional surfaces dual 



L. von Smekal / Nuclear Physics B Proceedings Supplement - preprint (2013) 1-42 



21 




SU(2) in 2+1 d 
with center vortex 



2d ising(A/xW) 
with interface 



-0 



Figure 15: Spatial center vortices share their universal behavior with 
interfaces in a spin model. Polyakov loops on either side of the vortex 
differ by the corresponding center phase. A thin vortex in SU(2) thus 
acts Uke an interface with antiferromagnetic bonds in the Ising model. 

to links at which the d.,, .5. in Eq. (50) are replaced by 
Ssi,sj+m mod q and are conveniently studied by introduc- 
ing analogous cyclically shifted boundary conditions, 

^it+giL = -Sf +Ci mod q , with c, = 0, 1, ...q-l . (52) 



Here, the interface free energies per temperature. 



Fj(^) = - In Rgit) 



-ln(Z,(a/Z,(0)) 



(53) 



are obtained from ratios Rq{c) of Potts model partition 
functions Zg(c) with cyclically shifted boundary con- 
ditions labeled by c over the periodic one, Z^(0). In- 
terfaces are suppressed in the low temperature ordered 
phase, and these ratios Rqic) tend to zero in the ther- 
modynamic limit. Interfaces in the spin model below 
Tc thus behave as the temporal center vortices in the 
ordered, deconfined phase of SU(A') above Tc- Com- 
plementary to that, in the disordered(confined) phase 
above(below) it is the interf ace( vortex) /ree energies 
that tend to zero such that the ratios Rq, R^ approach 1 
for all boundary conditions. Only at the critical temper- 
ature T = Tc do these ratios converge to non-trivial and 
universal values. 

In two dimensions these universal numbers. 



< < 1 , 



9.C 



= /?,(cO|j-^ , (54) 



for all cyclic boundary conditions c - (m,n) and all 
Potts models with 2"'^ order transitions, i.e., for q = 2, 
3 and 4, follow from the exact expressions of the cor- 
responding partition functions at criticality obtained in 
Ref. [105]. In terms of Jacobi elliptic theta functions, 

CO 

1?3(Z) = 2 ' 

«=-oo 

and with the definitions (for square lattices). 



TJm,n) = — —■^sixgin )&3iXgn ), where (56) 



Xq = l- 



s(V9/2) 



3 5 



:,1 , 



n [4' 6' 

for q = {2, 3, 4} , and // = T[l/4f/i2n^'^) , 



the partition functions of the ^-state Potts models with 

different cyclic boundary conditions at criticality can be 
expressed [105], for the Ising model with ^ = 2 as. 



Z<^°^= 472(4,4) -72(1,1), 
Z*':"' = Z^°:" = Zf + 272(2, 1) - 472(4, 2) , 



(57) 



"2,c 



"2,e 



2(1,1) ^ ^(0,0) _ 22(i'0) 

2,c 2,c 2,c 



for the q = 3 state Potts model as. 



Zf ^ = 673(6,6) - 373(3, 3) + 273(2,2) - 73(1, 1) , 



"3,e 

7(1,0) _ 7(0,1) _ 7(2,0) _ 7(0,2) 

^3,c ~ ^3,c ~ ^3,c ~ ^3,c 
7(0,0) 



(58) 



Z^J' -I- 3 73(3,1) -673(6, 2) 



7(1,1) _ 7(1,2) _ 7(2,1) _ 7(2,2) 

^3,c ~ ^3,c ~ ^3,c ~ ^3,c 

= zrV2-Z"''", 

3.C ' 3,c ' 



and the distinct ones for ^ = 4 as 



(59) 



7(1,0) 



: 274(2, 2) - 74(1, 1) -I- 474(4, 1) - 474(4, 2) , 



z^^f = 474(2,2) + 74(1, 1) - 474(2, 1) , 
^(m ^ _27,^(2, 2) - 74(1, 1) + 474(2, 1) , 

= ^S"' - 4^(4, 1) + 474(4, 2) , 

7(2,2) 

4,c 



: 674(2,2) + 374(1, 1) - 874(2, 1) 



{n,m) 



y(m.n) 
^q 



Generally, the partition functions obey Z^ 
and Z'^-"'"* = Z*"'"'. The ratios vt-^f = Z^^f^/zff 
readily follow. 

These universal ratios can be used to determine the 
critical couplings ySe of the transition in the 2+1 dimen- 
sional gauge theories with high precision [106], by re- 
quiring that fi - Pc sX Rk - Rq,c with q = N and the 
corresponding boundary conditions. For SU(2) in 2 + 1 
dimensions, for example, where, 

23/4 _ 1 

(60) 



R' 



(1,0) 

2,c 



1 



23/4 + 1 



and R: 



(1,1) _ 

2,c 23/4 + 1 



from the 2d square Ising model, this was used in [106] 
to determine the critical couplings of the deconfinement 
transition for lattices with up to A', = 9 links in the Eu- 
cUdean time direction with a precision typically two or- 
ders of magnitude better than previous Uterature values 
where they were available [99, 107, 108]. In particular, 
the high precision allows to reliably determine the sub- 
leading 1 /A/rCorrections to the linearly increasing be- 
havior of jSc with Nf near the continuum limit in 2 -1- 1 
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dimensions, 

PciNt) = 1.5028(21) A?, + 0.705(21) 

-0.718(49) 



(61) 



N, 



The slope of the leading term determines the critical 
temperature in units of the dimensionful continuum 
coupling ^3 of the (2+l)ii theory, yielding T^/g^ - 
0.3757(5). By standard arguments, from this result one 
can also read off the temperature dependence of the cou- 
pling at a given A^, [106], 

T, 0.270(2) t ^ , 
g^ Nt Ic I + 1 



t = TIT, - 1 



(62) 



For each fixed A',, we thus have precise control of the 
temperature T - \ /(aNt), where a is the lattice spacing, 
by varying the lattice coupling /3. The physical length 
of the spatial volume then follows from 



1 



N,T N,T, 1 + t(J3) 



(63) 



3.2. Twisted boundary conditions, electric and 
magnetic center flux 

As mentioned above, the different choices of twisted 
b.c.'s in pure SU(A^) gauge theories at finite tempera- 
ture are labeled by two Zftr-valued vectors, m and k, and 
fix the total numbers of vortices mod A^ through the or- 
thogonal planes of the Euclidean IjT x Lf spacetime 
box as sketched in Fig. 16. There are did - l)/2 purely 
spatial planes and the did - l)/2-dimensional vector 
m denotes the total conserved, Z^y-valued and gauge- 
invariant magnetic flux through those, as generated by 
a static center monopole. The ones through the d tem- 
poral planes labeled by the (/-dimensional vector k are 
the universal partners of the Potts interfaces. 

In order to understand the inequivalent choices for 
imposing boundary conditions on the gauge fields A, 
which are blind to the center of SU(A^), one first 
chooses Aix) to be periodic with the lengths of the 
system in each direction (1 up to gauge transformation 
Q^^."^) £ SU(A^), which is physically equivalent to peri- 
odic boundary conditions, 

A,ix+L^fi) = a""(x) = a^ix) (A,(x) - id,) i^^ix) . 

(64) 

In the 1/r xU' box at finite temperature, we use 





3d~line defect: Za-vortex, maps circle to a iion-contractible 
loop in RP(2), same happens in SO(3) = SC/(2)/Z2. 



Figure 16: Zo-vortex, e.g., as in nematic liquid crystals or pure 
SU(2) gauge theories (left). Center vortices in the + 1 dimensional 
gauge theory lower their free energy by spreading out: magnetic ones 
through the d(d - l)/2 spatial m-planes spread at all T, while those 
thi'ough the d temporal K-planes get squeezed with temperature (right). 

Then, compatibility of two successive translations in a 
(yu, v)-plane requires that (no summation of indices) 

Q.^ix+LyV)Q.yix) - Z^yQ.yix+L^^fl)Q.f^ix) (66) 
with Zf,y - e^m 'v. /A' ^ and n^v - -ny^, e ■ 

The total number modulo A^ of center vortices in a (/v, v)- 
plane is specified in each sector by the corresponding 
component of the twist tensor n^y. In d - 3 spatial 
dimensions magnetic center flux m through the box is 
given by n,-; = eijkmk, in li = 2 by n,j = e,j m. The time 
components no; = define the temporal twist k € Zf^. 

With these inequivalent choices of boundary condi- 
tions, the finite-volume theory in 3 -(- 1 dimensions de- 
composes into sectors of fractional Chern-Simons num- 
ber iv + k ■ m/N) [109] with states labelled by \k, m, v), 
where the integer v is the usual instanton winding num- 
ber These sectors are connected by homotopically non- 
trivial gauge transformations Q.[k, v]. 



Q[F, v'] |k, OT, v) — \k+k' ,m,v+v') . 



(67) 



Note that such a gauge transformation must be multival- 
ued in SU(A^) and hence singular to change the twist. 

A Fourier transform of the twist sectors Zkik,m,v), 
which generalizes that for 0-vacua as Bloch waves from 
v-vacua in two ways, by replacing v ^ iv + k ■ m/N) 
for fractional winding numbers, and with an additional 
Z?,-Fourier transform w.rt. the temporal twist k, yields. 



Zgie, fh, 6) 



Jcij{k,v) 



Zkik,m,v) . 



(68) 



k,v 



Up to a geometric phase. 



L^^=/3^l/T , and L,- = L, i^l,...d . 



(65) 



ioik, v) = Ine ■ k/N + Qiv+k ■ m/N) , (69) 
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the states in the new sectors are then invariant imder the 
non-trivial Q.[k, v] also, 

Q.[i v] 1^, M, 0) = exp{-ia>it v)} M, 6) . (70) 

Their partition functions Zg are classified, in addition to 
their magnetic flux iti and vacuum angle 0, by their Z^- 
valued gauge-invariant electric flux in the ^-direction 
[90, 110]. Here, we do not consider finite and omit the 
argument 6 in the partition functions henceforth. Mag- 
netic flux is irrelevant for the deconfinement transition 
also, because for aU combinations of magnetic flux fh 
with either electric flux ^ or temporal twist k [104], 

Zg(e, m) — > Ze(e, 0) , and Zi^(k, m) — > Zj^(k, 0) . 

(71) 

There is no analogue of magnetic center flux in the 3d 
spin systems, and center monopoles always condense 
because they are screened by the spatial string tension 
which is non-zero on either side of the deconfinement 
transition, likewise. Different combinations of temporal 
£and spatial m-twists, with difl'erent fractional topolog- 
ical charge, can however be used to measure the topo- 
logical susceptibihty without cooling [111]. 

The essential structure here is, however, that the tem- 
porally twisted vortex ensembles Zk(k) in the d + 1 
dimensional pure SU(A') gauge theory are related to 
those with fixed units of electric flux Z^i^ (and a per- 
haps more intuitive physical interpretation) by a d- 
dimensional discrete Z^r-Fourier transform [90], 

^^(^ = Y^Tj e^'^'^*>Z,(^) , (72) 

no matter what the magnetic flux is, which is why the 
argument w! is dropped here and in the following also. 

The role of the electric flux ensembles is best under- 
stood in terms of the translationally invariant flux be- 
tween a fundamental color charge at some point x in the 
finite volume and its mirror charge at j? + in a neigh- 
boring volume in the direction of the flux ^ [112], 

^ = -(trto/^'(;?+^L))) , (73) 

Z^(0) ^ V/no-flux' ^ 

where the Pi^'s are untraced fundamental Polyakov 

loops, including any potentially non-trivial transition 
function Qoix) in the time direction. With time-ordering 
from right to left, 

P(^ = Q.l(0,£)T[exp[i I Ao(t,^dt]). (74) 



The gauge field att =/3 difl'ers from that at ? = by the 
gauge transformation Qo(0, x) and we must undo this 
to make tr P{x) gauge invariant also under non-periodic 
gauge transformations g{x) which change the transition 
functions as 

o;(x) = gix + L^fi) 0^(x) g^x) . (75) 

With this, the Polyakov line in Eq. (74) transforms as 

P'i£)=giO,^P(^g\0,^, (76) 

and its trace is invariant under such gauge transforma- 
tions even when g(J3, x) + g(0, x). 

In order to derive Eq. (73) we first decompose the 
total electric flux vector ^ = ej^t, with e, - 1, ...N- 1, 
and perform successive translations by L in the direction 
of the individual unit vectors For each step we use 
the boundary conditions (64) for the gauge fields and 
the cocycle condition (66) for the transition functions 
Q.^, to show that in a fixed twist sector with boundary 
conditions k, and Z,o - e^^'"'^'^'^, 

P{x + Led = i^S(0, X + Lei) x) P{x) n;'"(0, x) 
= Z,o Q,(0, fi,')(0, f) P{^ O,!"(0, 
= Z,oQ,(0,f)P(x)Q,!"(0,f) (77) 

where we have used P{x) - T( expi A()(t,x)dt) to 
denote the normal Polyakov line as suitable for periodic 
boundary conditions. If we repeat this step for an arbi- 
trary translation by Le e iZ^, we pick up a product of 
center elements for the twists along the way, 

Pix + 11)= e'^"' ^-^1^ Q(0, P(f) (0, f) , (78) 

and Q(0, x) stands for the product of spatial transition 
functions when going from x to x -H Le. Changing the 
order of two spatial steps in different directions, L^i <-> 
Lej in the path from x to x + Le, multiplies Q. by the 
corresponding spatial twist Z, J. Because Z^Z^ = 1, their 
order does therefore not matter in Eq. (78). We have, 

Pix) Q'"(0, X) pHx + Le) Q(0, x) = e2;r;?-l/iV j ^ ^^g-, 

in an ensemble with fixed temporal twist k. When all 
boundary conditions are spatially periodic, Q(0,^ = 1. 
We can then insert the unity in Eq. (79) into the parti- 
tion functions with fixed boundary conditions and sum 
over all temporally twisted sectors to establish Eq. (73) 
from Eq. (72). The subscript 'no-flux' in Eq. (73) indi- 
cates that the expectation value is meant to be taken in 
the enlarged ensemble Z(,(0) = X]^Z^(^)/Af'' correspond- 
ing to this sum. If there is no magnetic flux, rh = Q, 
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we can always apply time-independent (non-periodic) 
gauge ti-ansformations with 

g(f + Ui) = gi^ ajio, f) ^ o;(o, f) = 1 . (so) 

We can thus assume without loss that Q(0, ^ = 1, as in 
Eq. (73) which then remains manifestly invariant under 
spatially periodic gauge transformations in this form. 
When i= we can still use Eq. (73) for a single di- 
rection, but for diagonal electric fluxes in a plane with 
magnetic twist we must include the combination of spa- 
tial transition functions Q(0, ^ in the expectation value 
on the right as in Eq. (79) which is then, however, invari- 
ant under non-periodic gauge transformations as well.* 
Eq. (73) emphasizes the physical interpretation of the 
electric flux sectors, as providing the necessary mirror 
charges for the sectors of SU(N) with non-vanishing to- 
tal color charge, dual to those with fixed boundary con- 
ditions. Unlike the standard product tr(P(x)P' (x + Ret)) 
of static charges at distance R in a periodic ensemble, 
which is not gauge invariant, the electric fluxes (73) de- 
termine the gauge-invariant color-singlet free energy or 
potential between static charges at distance L\e\. They 
contain no ultraviolet divergent perimeter terms and no 
short-distance Coulomb contributions either. 

5.5. Electric fluxes and selfduality ;n 2 + 1 dimensions 

Independent of universality and scahng, exact maps 
between the spin systems and their dual theories, in 
terms of disorder variables on the dual lattice, are pro- 
vided by Kramers- Wannier duahty [113]. In 2 dimen- 
sions, the dual of a spin model is a Zjy spin model 
again. In particular, the 2d (j'-state Potts models (50) are 
selfdual for all q. In 3 dimensions, the spin models are 
dual to Ziv gauge theories. The best known example is 
the 3c/ Ising model whose dual partner is the Z2 gauge 
theory, and this system is relevant for the dual universal 
behavior of center vortex and electric flux ensembles in 
the 3 + 1 dimensional SU(2) gauge theory [112, 114]. 
In 4 dimensions, the dual of a Z^r gauge model is a Zjv 
gauge model again [113]. 

One particular aspect of the SU(A^) gauge theories in 
2 -H 1 dimensions with N - 2, 3 and 4, that has been 
overlooked until recently, follows from the selfduality 
of the corresponding 2d (j'-state Potts models and uni- 
versality. The center vortex and electric flux ensembles 



*This slight subtlety was overlooked in [1 12]. To verify the invari- 
ance of Eq. (79) under general, non-periodic gauge transformations is 
a straightforward exercise in manipulating transition functions using 
their cocycle condition (66) and transformation law (75). 



of the 2 -I- 1 dimensional gauge theories are mirror im- 
ages of one another within the universal scaling window 
around criticality [100, 115, 116]. This is obvious from 
the fact that, as we will see, the exact duahty transfor- 
mation of the 2d g'-state Potts models in a finite vol- 
ume is precisely of the structure of 't Hooft's relation 
between the temporal center vortex and electric flux en- 
sembles in Eq. (72) discussed above. 

The duahty transformation and the selfduality of the 
^-state Potts models in 2 dimensions has long been 
known for infinite systems. A particularly simple proof 
based on the random bond-cluster representation was 
given in [93]. Far less is known in a finite volume with 
translationally invariant boundary conditions, however, 
where ensembles with different boundary conditions in 
general mix under duality transformations. One case 
where this is known is the 3d Ising/Za-gauge theory sys- 
tem [117-119]. For the Potts model, using the random 
cluster methods developed in [105], we were able to ob- 
tain the following exact duality transformations for all 
^ on a finite 2d torus as discrete 2d Fourier transforms 
over all ensembles Z^™'"' with cychcally shifted bound- 
ary conditions ? = (m,n) [115], 

' ^ m,n 

m,n,r,s = 0,1,.. -q-l , (81) 

where A^sites is the total number of sites of the 2d lattice 

(i.e., A^sites - onanN X N squarejattice), K - J/T 
is the coupling per temperature, and K its dual obtained 
from 

(e" - l)(e^ - 1) = ^ , (82) 
as usual, with temperature mirrored around criticality at 

K = K = Kc= ln(l -I- V^) . (83) 

For q - 2, and with K K/2, conventionally, this du- 
ality relation reduces to an analogous known result for 
the Ising model [120]. With q = 3 it agrees with a re- 
sult obtained for the planar or vector Potts model 1 121] 
which is equivalent to the standard one in that case. 

Before we turn to the sketch of the proof, note that 
the structure of the finite-volume duality ttansformation 
in the Potts models (81) is precisely the same as that 
in the relation between the temporal center vortex and 
electric flux ensembles in Eq. (72). The temporal center 
vortices are the universal partners of the interfaces in the 
spin model, c.f., Fig. 15, and the electiic fluxes are their 
duals, with d = 2 again universally related to ensembles 
with interfaces, but at the dual temperature, swapped 
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around criticality. In J = 3 spatial dimensions the pat- 
tern is the same [1 14], the center vortex free energies in 
SU(2) share their universal behavior with those of the 
interfaces in the 3d Ising model, and electric fluxes cor- 
respond to ensembles of the dual 3d Z2 gauge theory 
with anti-periodic boundary conditions. But they are 
not exact mirror images of one another in <i = 3. 

The general formula (81), for all <7-state Potts models 
on a 2d torus, was first given to my knowledge in [1 15]. 

I have not pubUshed its proof nor have I seen it in the 
literature before. So I will sketch it here, in the rest of 
this subsection which will not be needed later on and 
which the reader may thus choose to skip ad libitum. 

The proof uses mostly standard techniques [93, 122], 
and it is based on the setup and results of Ref. [105]. 

We start from the Hamiltonian (50) without external 
field {H = 0) but with cyclic interfaces such that 

-p^M =kY, S(Si, G^;^'"^sj) , (84) 

where the notation is such that we have introduced ma- 
trices G*^™'"' on the links connecting nearest neighbors 
and 6{s,s') = 6s,s' = 1 when s = s' and zero oth- 
erwise. For the cycUc interfaces we introduce straight 
lines (seams) Lj and Ly in the x and y directions on the 
dual lattice and use 

r 5-, {ij)€L;, 

G^;'"> = I S", {iJ}€L;, (85) 
{ 1 , otherwise , 

where S isaqxq matrix for an elementary cyclic shift. 
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such that along the x-direction the spins get shifted m 
times at the links dual to {i.e., crossing) Ly and n times 
in the y-direction at those dual to L^. The cyclically 
shifted boundary conditions implemented in this way 
are special cases of the general permutations of spin 
states across the boundary considered in Ref. [105]. We 
will use the random-cluster representation of the Potts 
models on a finite torus from there. With 

V = - 1 (87) 

the random-cluster representation of the partition func- 



tion with (m, n)-boimdary conditions then follows from 

Z^-«> = exp {k J] S(s^, G^J'-'sj)} (88) 

i.v) m 

is] (ij) graphs 

This is the standard random-cluster representation of the 
infinite Potts models [93, 122, 123]. It is vaUd also on 
the 2d torus in this form as long as we only consider 
the cychc shifts of the spin states across the seams in 
Eq. (85) as will be shown below. The product in the sec- 
ond line is multiplied out and the spin sum represented 
as a sum over all graphs where a hnk with aligned spins, 
6(Si,G\j'"^ Sj) = 1, is said to be occupied while those 

with 6(si,G^^'"\j) = are empty. The possible con- 
figurations are then identified with graphs consisting of 
clusters connected by occupied links, whereby each dis- 
connected site counts as its own cluster. So all sites are 
part of a cluster. Each occupied link (or bond) in a graph 
contributes a factor v, and their total number in a given 
graph is called Ni,. Each of the Nc clusters in a graph 
can be flipped independently as a whole and thus gives 
a combinatoric weight q to the particular graph in the 
partition sum. 

On the 2d torus one needs to classify the clusters 
by the possible winding numbers (cji.cji) of the non- 
selfintersecting closed loops around the torus that one 
can draw on them: 

- If one cluster extends across both boundaries in 
such a way that independent (0, 1) and (1,0) loops are 
possible, the graph is said to have a torus cluster (TC) 
of which it can have at most one. 

- A cylinder cluster (CC) extends across the bound- 
aries in such a way that only non-trivial loops with 
winding numbers (^1,^2) are possible, where toi and 
0)2 are coprime. We call the number of cylinder clusters 
in a graph A^cc. and they must all have the same pair of 
winding numbers. 

- And finally, island-type clusters are the ones on 
which only (0,0) loops are possible. A graph with 
only island-type clusters is called a torus lake (TL) be- 
cause the empty links form stacks which extend across 
the boundaries in both directions such that independent 
(0, 1) and (1,0) loops are possible on the dual lattice, 
perpendicular to these stacks. 

GeneraUy, the empty links on the original lattice X 
are dual to links which can be connected to form a clus- 
ter on the dual lattice X*. Therefore, a lake on the orig- 
inal lattice is a cluster on the dual lattice [105]. 

With this classification, and more general seams cor- 
responding to arbitrary and independent permutations 
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Gx and Gy of the spin states across each boundary in x 



The first sum is over all graphs with cylinder clusters, 
the second over all graphs with one torus cluster and the 
third over all graphs with only island-type clusters in a 
torus lake. See Ref. [105] for more details. 

In the second line (90) we have used that Gx - S "' 
and Gy = S" here, and that the weight of an {0)1,(02) 
cy Under cluster from [105] then reduces to 

CiGx,Gy;a>u(02) = tr(G^>G^^) = ti:[S'^'S"^) , 

(91) 

whereby, with the Kronecker delta modulo q, 

tr(5""*"5""*'^) = qdgimoji + noj2) (92) 




q , mcL)\ + na>2 = mod q , 
, otherwise . 



Moreover, unlike the more general case with arbitrary 
permutations Gi, Gy in [105], for the cyclic ones in 
Eq. (85) considered here, a torus cluster is only possible 
when m - n - 0, i.e., with periodic boundary condi- 
tions. But then it can have q different states just as any 
other island-type cluster, therefore 

( q , m = n = , 
T(m,n) = T(S'",S") = \ (93) 
I, , otherwise . 

This means that when we restrict to cyclically shifted 
boundary conditions, all clusters that are possible for 
such given (m,n) boundary conditions have the same 
weight, v'^^'^^s and this establishes that the infinite vol- 
ume form of the cluster representation in Eq. (88) re- 
mains formally valid. One only needs to figure out the 
allowed cylinder and torus cluster and torus lake graphs 
as in Eq. (90) which is the starting point for the duality 
transformation in the finite volume. 

The first task in obtaining the dual cluster represen- 
tation is to re-express the number of occupied Unks Nb 



and y directions the random cluster formulation of the 
Potts models on the torus as derived in [105] becomes, 



(89) 



(90) 



and the number of clusters Nc in Eq. (90) in terms of the 
corresponding numbers on the dual lattice X*. 

To achieve this, recall that for an arbitrary 2d lat- 
tice, the Euler formula relates the number of vertices 
V, edges E, faces F, and open ends A'open to the genus g 
of the surface on which it is drawn via the Euler charac- 
teristic, 

X = V-E + F = 2-2g-Nopen. (94) 

The Euler characteristic of the torus with g = I and 
without open ends, A^open = 0, is x = 0. Moreover, 
y = A'j is the number of sites of X, and E is the total 
number of links, that is the sum of occupied ones Nt 
and empty ones which is identical to the number A^^ of 
occupied ones on the dual lattice Thus E = Nh+Ny 
Finally, the number of faces is the number of sites on 
£\ i.e., F = n;. Therefore 

Nb=N, + N:-N;. (95) 

It also follows from the Euler relation on the torus with 
X = that the total number of sites Ns of the lattice 
X is given by the number of clusters Nc in a graph (re- 
call that every site is part of a cluster) plus the number 
of occupied hnks A^^, minus the number of independent 
non-homotopic (and non-selfintersecting) loops Ni on 
the clusters of the graph. Therefore, the number of clus- 
ters A^c in a graph is expressed as 

Nc=Ns-Nb+Ni. (96) 

All TC graphs, with one torus cluster on the original 
lattice, are graphs with one torus lake on the dual lattice 
and vice versa, i.e., (TC) = (TL*) and (TL) = {TC*}. 

On a TL graph, where all clusters are islands whose 
coasthnes have (0,0) winding numbers, one can walk 
around every pond which is a lake with (0, 0) coasthne 
inside a cluster. Therefore Ni = A'ponds- On a TC graph. 



{CO \ Q / (TL) 



Here this simplifies to. 



{CO 



{TO {TL) 
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one can walk around all ponds individually, but walking 
around all of them is homotopic to a point. In addition 
there are the (0, 1) and (1,0) loops, thus Ni - Appends + 1- 
Together with the fact that the islands of a TL graph 
are the ponds of its dual graph TC* and vice versa, this 
implies for 

{TL) = {TC*) : Ni - A?po„ds = A^,* - 1 , (97) 
{TO = {TL*) : Ni = ^ponds + \^K + \. 

The ponds become islands on the dual graph. So Appends 
is equal to the number of dual clusters A^* on TL* 
graphs, while A^* is the number of islands (and thus 
Appends) plus one on TC* graphs. 

If we insert the identities (97) together with the Euler 
relation (95) into Eq. (96), this leads to 

{TL) = {TC*) : Nc = -A^,* + N'l + N*, - \ , (98) 
{TC} = {TL*} : Nc = -N';+Nl+Nl + l . (99) 



CC graphs are selfdual. They contain an equal num- 
ber of cylinder clusters and cylinder lakes, and all of 
them have the same winding numbers (wi,W2). The 
dual graph is the one with clusters and lakes swapped, 
and it thus has the same number of cylinder clusters, 
A'cc = A'cc ^® same set of winding numbers. On 
a CC graph one can walk non-homotopically around ev- 
ery pond and along the (wi.wi) loop of every cylinder 
cluster, therefore, 

{CC} = {CC'} : 

A^; = Appends + A^CC = A^* lands + ^CC = K - 

Nc = -n: + n;+n: . (100) 

Now we use the relations in Eqs. (95), (98), (99), and 
(100), to re-express Nh and Nc in Eq. (90) in terms of 
the corresponding numbers A'^ and A'* on the dual lattice 
£,*. Labeling the graphs in the sums by their respective 
dual cluster configurations, this yields for the partition 
function with (m, n) boundary conditions. 



Y, - q'^'d.ijnu, +na>2) + J„,o^„,o J] ^ q'^^' + ^ ^ , (101) 

ICC-) {TL'l {TO > 

and its discrete Fourier transform over m,n = 0, ... ^ - 1 becomes, 

{/ \''^^ 
^iKiirm^snVg m qK-^ Sq{mu, + na>2) (102) 
...... m,n {CC*) 

- Z (!) v., (!)">-= }^ 

{TL*) \ ' {TC*) ^ ' ' 



We have factorized another q in Eq. (102) to get the right 
power A^* in the middle term for which the sums over m 
and n collapse to the single m - n - Q term. The ad- 
ditional factor q^ in front of the last term in Eq. (102) 
arises because the TL graphs are independent of the 
boundary conditions m and n and therefore their sums 
yield qSy^ and qSsfi^ respectively. In the first term, the 
graph sum extends over all kinds of graphs with cylin- 
der clusters for all pairs of coprime windings a>\ and ai2, 
and the periodic Kronecker delta. 



q-\ 



^ (=0 

selects the right ones for each pair of boundary condi- 



tions im,n). The duahty transformation is thus com- 
pleted upon realizing that for any pair of coprime 
{0)1,(02), we can write. 



2 



(104) 



m, n t=Q 



= ^ qdqiojyt - r) q6q{p}2t - s) = q^6qUo2r - ojis) . 

The last step might not be immediately obvious. It is 
verified by showing the following:' 

thank Nils Strodthoflf for verifying that tliese two conditions are 
indeed satisfied. 
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(a) There is at most one ^ in {0, ... ^ - 1} which solves 
the two conditions cjit - r mod q and (D2t = s mod q at 
the same time. If it does, then the condition W2?"- wis = 
mod q is also satisfied. 

(b) If there is no solution t to the first two conditions, 
then the new condition a)2r - ojis - mod q does not 
have one either. 

Therefore, we have shown that the terms in brackets 
on the right in Eq. (102) agree with the random cluster 
representation (90) of a Potts model partition function 
Zg ' with boundary conditions (-s, r), and with the re- 
placement V q/v which is equivalent to replacing the 
coupling per temperature K by its dual K according to 
their duality relation in Eq. (82). Analogously rewriting 
the prefactor on the right in Eq. (102) with Eq. (82), 

v^:=(.^-I)'v; = ^!L_, (105) 
(e^ - 1)^' 

and rearranging Eq. (102) to solve for the terms in 
brackets on the right, equalling ' (K), then finally 
yields the selfduality relation (81) of the 2d ^-state Potts 
models (50) for all q. The result here is actually a bit 
more general than that given in (8 1) where we have used 
A^sites = Ns = N* for a 2d square lattice. More generally, 
the duaUty relation valid for any lattice X on a 2<i torus 
without open ends, and symmetric in X and £,*, here 
follows from Eq. (102) as 

(e^ - 1)"'^' Z'^-'-'-\K) = (106) 

_ ]^y!^s 1 ^ ^lni(rm+sn)lq Z^"''"\K) 

At criticahty, with K = K, \hc prefactors cancel and 
one verifies that the values of the partition functions 
for the various cyclic boundary conditions in Eqs. (57), 
(58), and (59) reproduce themselves under the discrete 
Fourier transform (106) for q = 2,3 and 4 as they must. 
The finite-volume duality relation is valid, however, for 
all q, including q> 4 with first order transition. 

In the 2-1-1 dimensional gauge theories, temperature 
is the same on both sides of the ZA?-Fourier transform in 
Eq. (72). As a consequence of the selfduality of the cor- 
responding spin models, however, the free energies of 
spatial center vortices and those of the confining elec- 
tric fluxes are mirror images of one another within the 
universal scaling window around the second order phase 
transition. 

3.4. Universality and finite-size scaling 

The ratios of partition functions with (^, OT)-twisted 
and periodic boundary conditions in the SU(A^) gauge 



theories are obtained from multiplying coclosed stacks 

Q.*(fi,v) of plaquettes by the center elements such 
that one plaquette in every plane of orientation (jj., v) is 
multipUed by the Z^y = e^"'"""^^ corresponding to the 
twist tensor n^y (c./, Eqs. (66)), 

Z^ilA) ^ / UdU exp {-PY.^ 1 Re(Z^ tr t/J} 
ZkM /n^/f/exp{-j82jRe(trf/^)} 

^uA \ u (107) 
I, 1 , otherwise . 

These coclosed stacks of plaquettes 0*(yU, v) are dual to 
non-homotopic closed surfaces or lines on the dual lat- 
tices in3-i-l or 2-1-1 dimensions, respectively. 

If the total number of plaquettes in all stacks Q*(/U, v) 
is A'^, the ratio of partition functions in (107) can be 
converted into a product of A'^ ratios, each of which rep- 
resents the expectation value of the «* twisted plaquette 
2W[/(n) in an ensemble Z^""'^ with n-l already twisted 
ones, where n runs lexicographically through the pla- 
quettes of one 0*(yu, v) stack after another, to implement 
twisted b.c.'s in all (/i, y)-planes according to the non- 
zero components of the twist tensor n^y, 

^^(^'"^) = n ^ = n (eX^'-^'x% 

Zk(0,0) ^ i 7(«-l) i i \ An-D' 

n=l n=l 

zf = z,(o,o), zl""^ = Zk(la) . (108) 

In this way, one ratio of partition functions requires A^^ 
independent Monte-Carlo simulations [112, 124]. 

For the temporal fc-twists one introduces coclosed 
stacks Q* aligned with the time direction between one 
pair of adjacent time-slices, where they are dual to spa- 
tial lines in 2 -H 1 and surfaces in 3 -H 1 dimensions, just 
as the respective spin interfaces in 2 and 3 dimensions 
as illustrated in Fig. 15 above. 

Once the reduced temperature t(fi) is determined in 
terms of the lattice coupling /3, as from Eq. (62) for 2 + 1 
dimensional SU(2), and with the spatial lattice size L 
from Eq. (63), a finite-size-scaling analysis can be per- 
formed. In the vicinity of the 2"^* order phase transition 
generalized couplings such as the vortex-ensemble ra- 
tios Rk(k), for sufficiently large L, only depend on the 
ratio of L and the large correlation length ^ which in the 
infinite volume diverges as 

^ = /i?r + ■■■ , f ^0* . (109) 

with the correlation-length critical exponent y, where 
for the 2d Potts models v = 4/3 for ^ = 1 (percola- 
tion), V = 1 for ^ = 2 (Ising model), v = 5/6for^ = 3 
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Fj{N,K) = -\nRl'll\N,K) 




-10 -8 -6 -4-2 2 4 



x = Nt 

Figure 17: The universal scaling function fi(x) calculated for a (1,0) 
interface corresponding to one antiperiodic direction from the ratio 
exact finite volume partition functions Z^LiVz^'i'j' of the NxN square 
Ising model in Ref. [125]. 

and V - 2/3 for q = 4 [122] . Here we are particularly 
interested in the ratios 7?^"' "' which, up to finite size cor- 
rections ~ A^"" on an X square lattice, depend on 
X - N^l^t in a universal way. The corresponding univer- 
sal scaling functions fi{x) for the interface free energies 
per temperature, suppressing the indices (m, n) for the 
boundary conditions, are given by 

F,(N,K) = MN^'^t) + ciN'" + ■■■ , with 
//(O) = -In/?, , and (110) 
fi(x) = o-o (-xf + ■■■ , x^ -oo , 

where K - KJ(t+l), Kc = ln(l -H yfq), and Rc stands for 
the universal ratio /?^'".'"' with the particular combination 
of boundary conditions that is being used in the q - 2, 
3 and 4 state Potts models as per Eqs. (57), (58) or (59). 

For the 2d square Ising model, with v = 1, the fi- 
nite volume partition functions Z^"^~\n,K) have been 
obtained exactly for all combinations of periodic and 
anti-periodic boundary conditions in Ref. [125]. We can 
use these results for the ratios R^^f^ ™d R^^^2 to obtain 
the corresponding universal scaling functions //(x) from 
Eq. (110). In Fig. 17, we plot F,{N, K) for the (1 , 0) in- 
terface, for example, over x - Nt for = 100 and 
A^ = 1000. We have tested that the diff'erence between 
F,{N,K) and f,{x) vanishes as A^"'^ with w = 2 [106]. 
The corresponding small residual finite-volume efifects 
are observed, e.g., in the deviation of the A^ = 100 result 
around x - -5 from that for A^ = 1000. To check which 
lattice size is needed for any desired accuracy, one can 
calculate the deviations of Fi{N, K) from the asymptotic 
slope (To for x — > -oo which is also known analytically. 

Generally, the interface tensions (per temperature) in 



the d dimensional spin models are defined in the ordered 
phase as 

a-i^ \imN'^'''^'^F,{N,K) , K > K, . (Ill) 

N—tco 

Comparison of this definition with Eqs. (110) thus en- 
tails that, if CT/ is finite in this limit, the interface ten- 
sions near criticality must behave as 

cr, - o-Qi-tf + ■ ■ ■ , with /J - (d- l)v , t ^ 0' . 

(112) 

The exponent of the interface tension // is tied to the 
correlation length exponent v by one of the so-called 
hyperscaling relations. In particular, for the 2d square 
Ising model with one antiperiodic direction, with fj. - 
V = 1 and thus ctj - -o-qI + ■ ■ ■ , the asymptotic slope 
(To follows from Onsager's famous result [126], valid 
for all < r < in the ordered phase with A^ — > 
oo, for the tension of a straight (1,0) interface in the 
thermodynamic limit, 

^(1,0) ^2/:-Hlntanh/:. (113) 

When the spin-coupling per temperature K - J/T is 
expanded about criticality at Kc - ln(l H- V2)/2 this 
leads to,^ 

cr[,''°^ = 21n(l + V2) , (114) 

and this agrees with the asymptotic slope of fi(x) for 
X — » -oo in Fig. 17 as obtained from the exact finite 
volume partition functions 2^^;°' ^^d Z^^;^'' of Ref. [125]. 

Many alternative and simplified derivations of On- 
sager's result (113) were given since then and can be 
found in text books [122, 127]. A particularly interest- 
ing one for us is given in Baxter's book [122], as it pro- 
vides the tension of a diagonal (1,1) interface. In our 
notations his result can be written as 

cr*'''> =21nsinh(2/:) . (115) 

As the free energies so far, these are all interface ten- 
sions per temperature. The proper free energy of an in- 
terface is given by -Tin 7?^, - TcrjN. The correspond- 
ing proper interface tensions, 

'L{e)= Oo-i , e^TIJ ^K-\ (116) 

per coupling / over temperature, for 9 - Q lo 6^ - K^^ , 
following from Eqs. (113) and (115) are compared in 
Fig. 18. At r = 0, every antiferromagnetic bond- 
coupling costs 2J in energy. The diagonal or zigzag 



We use the typical Ising model conventions for J such that K — » 
KI2 here as compared to the general Potts model formulae above. 
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S(6l) = 9ai 



Baxter (diagonal) 




kT 

T 



ln(l + V^) 



Figure 18: The exact interface tensions per coupling, 1.(0) = 6(Ti over 
temperature 6 = TjJ, for the straight (1,0) interface from Onsager, 
Eq. (1 13), and for the diagonal (1,1) interface from Baxter, Eq. (1 15). 

interface is twice as long as the straight one. It needs 
twice as many antiferromagnetic couplings, hence, in 
units of y, E(0) = 2 for the straight and 4 for the diag- 
onal interface. In the isotropic limit, when the correla- 
tion length is much larger than the lattice spacing, on the 
other hand, the underlying lattice does not matter any- 
more and the diagonal interface through a large finite 
volume is only a factor of V2 longer than the straight 
one, i.e., from Eq. (115) around K - Kc, 



2V21n(l -H V2) = V2c 



.(1,0) 



(117) 



The same square-root ratios near are observed for the 
3(i Ising interfaces [128], 



(1,0,0) (1,1,0) 



Ji,i,i) 





1 : V2 : V3 . (118) 



The spin-interface tension cr/ corresponds to the dual 
string tension & for spatial center vortices in 2 -i- 1 and 
spatial 't Hooft loops in 3 H- 1 dimensions for which a 
dual area law holds in the high-temperature Z^^-broken 
phase. The same square root ratios (118) are observed 
for the dual string tension in SU(2) [112, 1 14]. From the 
duality Eq. (72), on can show that they then must also 
hold for the electric fluxes in the confined ZAr-disordered 
phase below T^, and they are interpreted as the smoking 
gun of string formation: The free energy of orthogonal 
fluxes is minimized when their length is, i.e., when di- 
agonal strings through the volume form. 

At least for the 2+1 dimensional SU(2) gauge the- 
ory string formation just below the deconfinement tran- 
sition temperature is thus proven from universality and 
the observed square-root ratios of the exact Ising model 
interface tensions in the isotropic limit. 

Unlike the spin models, however, whose interface 
tensions at T = follow ratios 1 : 2 : 3 as for isotropic 
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Figure 19: Vortex free energies Fk in SU(2) above for Nt = 4. 



flux, the electric fluxes in the SU(A^) gauge theories are 
expected to also show the square-root ratios for string 
formation at zero temperature. The numerical evidence 
from early Monte-Carlo results on the lattice is aflirma- 
tive that this is indeed the case [129]. 

The {2+\)d SU(2) vortex-ensemble ratios Rk(k) were 
calculated on N,y.N^ lattices for A^, = 2 to 10, each with 
spatial sizes up to Ns - 96 and in a suitable window of 
lattice couplings (3 around fi^ to test the finite-size scal- 
ing (ESS) in Refs. [106, 115, 116]. Eor each fixed N, 
one generally observes very good scaling of the avail- 
able data for all Ns when plotted over the ESS variable 
X - L^^^t. With V = 1, the physical spatial length L 
from Eq. (63) in units of T^', and the reduced temper- 
ature t(J3) from the lattice coupling p via Eq. (62), we 
define. 



X = LTa 



N 1 + m 



(119) 



Eor the resulting vortex free energies F^, for k = (1,0) 
and (1,1), we then obtain accurate one-parameter fits 
via 

Fk(x) = M-Ax) , (120) 

to the exact universal scaling functions //(x) for (1,0) 
and (1, 1) b.c.'s computed from the results of Ref [125] 
as in Eig. 17. This determines the single non-universal 
parameter A = A(N,) which is dimensionless and relates 
the SU(2) ESS variable to the Ising one. 



J^Ising - - A(N,) Xsi](2) 



(121) 



whereby the minus sign reflects the interchange of high 
and low temperature phases between the two. 

How hard it is to precisely determine the dual string 
tension cr directly from the data for the vortex free en- 
ergy is shown in Eig. 19. We need to get to asymptoti- 
cally large x, c.f., Eqs. (1 10), but that requires larger and 
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Figure 20: The ratio Ri,(k) of tlie partition function with twist k = (1,0) over the periodic ensemble in the 2+1 dimensional SU(2) gauge theoiy 
compared to Re(^ for one unit of electric flux = (1,0) relative to the no-flux ensemble (left), and to its mirror image (right). The data here is 
obtained for N, = A and spatial volumes up to Ns = 96 [115]. The dashed lines represent the universal scaling function from the 2d Ising model, 
exp{-//(-/ix)) for Rt and exp{-/)(/lA-)) for from selfduality. 



larger L - aNs in order to stay within the universal scal- 
ing window for sufficiently small t. Note that x - 4 here 
roughly amounts to x » -6 in Fig. 17 where A^^ - 100 is 
not large enough to suppress finite-size corrections even 
in the Ising model. The computational costs in the d+l 
dimensional gauge theory roughly increase as 
so it becomes rather expensive to beat the FSS correc- 
tions to the dual string tension by brute force, and this 
was observed even more so in = 3 [112]. 

Here, our one-parameter fits to the vortex free ener- 
gies as per - fj{-Ax) to the exact universal scal- 
ing functions help tremendously because we can ob- 
tain accurate data at comparatively low computational 
costs for these fits from small x values, i.e., small vol- 
umes. Once A is fixed, however, we can compute the 
dual string tension and, as we will see, also the string 
tension for electric flux from these. 

This is because with the hyperscaling relations /i = 
(d - l)v, the product of interface tension and correlation 
length, cTi^^'^^ is independent of t near criticality. In 
fact, this product is another example of a universal am- 
plitude ratio [8]. To be precise, if we use the so-called 
exponential correlation length in the disordered phase 
of the spin model above T^, ^g^p, then the universal con- 
stant is given by 



where 



(122) 



r: 



1 , for q - 2, 3, 4 , in d - 2 , 
■■^gap - I 0.40(1) , for ^ = 2 , in = 3 . 

The d = 2 results are again exact [93, 122, 127],'' while 

'The q = 2 result is found in any of the textbooks, for = 3 it 



the 3c/ Ising ratio is determined numerically, e.g. , see the 
review in Ref. [8]. 

The exponential correlation length corresponding to 
^gap on the gauge theory side is that of the Polyakov 
loops in the confined phase below T^, which is given by 
the string tension cr and temperature as 



1 

r 



cr 
T 



(123) 



Therefore, we have a universal relation between the 
string tension per T below T^, playing the role of ^g^p ' 
in Eq. (122), and the dual string tension & above Tc, as 
the interface tension cr/ in Eq. (122), 



(T = Rta^pi^/T^y' around criticality, or 

a- = pTli-ty + ■■■ , T <Tc, and (124) 



-=<gap 



d-lJd-l)v 



,T>T,, 



where we have introduced an unknown non-universal 
constant p which drops out from the universal ratio, and 
we have replaced one factor T by Tc in the string tension 
in second line as we work at leading order in the reduced 
temperature t here. Instead of determining the constant 
p from the vortex free energies in the high temperature 
phase at asymptotically large x - {TcD}^^ t, here we 
can use the exact universal scaling function //(jc) from 
the Ising model: We know the asymptotic slope of //(x) 
which is 2 ln(l -i- V2), therefore, from 



Fi,(x) = fi{-Ax) p = 2/lln(l -H V2) 



(125) 



follows from the hard hexagon model [122], and for q = A probably 
from the Baxter-Wu (triangular three-spin) model [122, 130], but I am 
unsure about the status of a mathematical proof of the latter. 
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variable Ax for N, up to 10, and N, up to 96 as in Fig. 20, from [116]. pigu^g 22: Selfduality in 2 + 1 dimensional SU(3) from [100]. 



Once we extract A from our one-parameter fits at small 
X where we have very accurate data, this then equally 
accurately determines the string tension and its dual 
around Tc- Moreover, we can determine this single non- 
universal parameter AfoiN, = 4, 5 ... 10 which allows 
a polynomial fit, 

A(N,) = ioo + b/N, + c/Nf , (126) 

with the extrapolated result = 1.354(25) [115]. This 
then determines the leading behavior of the continuum 
string tension and its dual around the phase transition, 

o- = Ao,T^2 ln(l + V2) |f| + ■ ■ ■ , f 0" , and 
<5- = ioor,21n(l -H V2)f+--- , f^0+. (127) 

Once we have fixed the single non-universal parame- 
ter A we can plot the calculated center vortex free ener- 
gies Fk or the corresponding ratios Ridk) in comparison 
with exp{-//(-/lx)) as in Fig. 20 for k = (1,0). The data 
for all different spatial lattice sizes collapse beautifully 
onto the universal curve, the finite-size scaling descrip- 
tion works rather well. Together with the corresponding 
k - (1,1) result for two orthogonal center vortices we 
can perform the 2d Z2-Fourier transform (72) to obtain 
the ratios of electric flux e over the no-flux ensemble in 
Eq. (73), 

^''^^"^ ^IM^^ (tr(P(f)pt(f + e-L))l^ ^^^ . (128) 

These are equally well described by exp{-//(/lx)) with- 
out refitting the non-universal parameter A, they are mir- 
ror images of the Ridk) in a surprisingly large scaling 
window around criticality as seen in the right panel of 
Eq. (73), where we swapped jc — > -jc for the electric 



fluxes to demonstrate that. This is the manifestation of 
the selfduality of the corresponding 2d spin model. 

Moreover, if we rescale the ESS variables by the pa- 
rameters A{N,) determined from the one-parameter fits 
for the individual A^, = 4 , ... 10 data sets, the scaling 
of the data over all spatial lattice sizes A^, can be ex- 
tended to include different A^, lattices. With this renor- 
malization, one thus essentially obtains good contin- 
uum results akeady with the coarsest N, - 4 lattice 
here, as seen in the zoom-in plot around criticality in 
Fig. 21. The rescaled x-range there, roughly corre- 
sponds to X = [-1.33, 1.33] in Fig. 20, and the dashed 
line is obtained from the universal scaling function as 
exp{-//(x)) and thus now independent of A^,. The slight 
deviations of the vortex free energy from this curve near 
X = I seen here, are predominantly due to the leading 
finite-size corrections ~ L analogous to what is be- 
ing observed in the ordered phase of the Ising model in 
Fig. 17 at negative x. 

The selfduality of the underlying Potts model is also 
observed in the 2d SU(3) gauge theory [100], as shown 
in Fig. 22. The dashed line in this figure marks the uni- 
versal ratio R^^'^^ from Eqs. (58). The two independent 
ones that one can form on a symmetric lattice for q = 3 
evaluate to 

= Rff = 0.30499982 . . . , and (129) 
M'"'^ = M^'^' = r[^'^^ = 0.19500018 .... 

i,c 3,c i,c 

Even though the data in Fig. 22 was obtained on a quite 
small N, -2 and A^, = 24 lattice, the intersection point 
of the vortex and electric flux ensembles Rk and Rg is 
almost right on the universal line. To exploit selfduality, 
where possible, turns out to be the fastest converging 
method to determine critical couplings by far 
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Figure 23: Critical couplings for (2+l)rf 5f7(2) (N, = 4) from self- 
duality compared to those of Ref. [106], with Ns up to 96 and the 
infinite-volume extrapolated result yS^ = 6.53661(13) shown as the 
narrow grey band here, from [116]. 

There is a long history of methods to extract criti- 
cal couplings or temperatures from simulations in fi- 
nite volumes, going back to using pairwise intersec- 
tions of Binder cummulants on successively larger lat- 
tices [131]. Hasenbusch later demonstrated that the ra- 
tios of partition functions with different boundary con- 
ditions could be used in the same way to obtain a much 
more rapid convergence with very good estimates al- 
ready from rather small lattices [132]. At criticality, 
these ratios tend to the universal values < Rc < I 
in the thermodynamic limit. In [106] it was therefore 
shown how to obtain critical couplings for gauge theo- 
ries from intersecting the ratios of finite volume par- 
tition functions with these universal fixed points, once 
their values are known. For {2+l)d SU(2) this led to an 
even much faster convergence than their pairwise inter- 
sections. At the time we thought this is the best method, 
but now look at the comparison of the results from [106] 
with the intersection points from selfduality in Fig. 23, 
again for SU(2) in 2 -H 1 dimensions as a benchmark. 
Already with A^^ = 16, and without any extrapolation, 
the result from selfduality is within the errors of the ex- 
trapolated best infinite volume result, jSc - 6.53661(13) 
from Ref. [106], as indicated by the narrow grey band 
in Fig. 23. 

The reason for this impressive result is that the lead- 
ing finite-size corrections do not change the position of 
the intersection point in They only move it upwards 
along a straight vertical line with increasing L in plots 
such as the one of Fig. 22. 

To see this explicitly, first use a finite-size scaling 
ansatz for the vortex ensemble ratios R/; as functions of 



the lattice coupling /? around criticality of the form 

Rk(J3) ^Rc + b{Ji- Pc) N]!' + cN-'^ + ■ ■ ■ , (130) 

Defining pseudocritical couplings j3c{Ns,Nt) in a finite 
volume by requiring that RkiP) — Rc for the method of 
[106] leads to, 

lic{Ns,N,) = li,{N,) - (c/b) A?;(-+i/>') + . . . . (131) 

These extrapolate to PciNt) from large spatial lattice 
sizes A^j at fixed numbers of time slices A^,. As a byprod- 
uct this method gives numerical estimates of the correc- 
tion to scaling exponent lo. 

With selfduality, however, one must then have 
Reifi) - RkiP) for like values of electric flux e and and 
temporal twist k at /3 = /5c- In fact, it is straightforward 
to verify that then, 

ReiJ3) ^Rc-b(/3- AO Nl'' + cN;'^ + ■ ■ ■ , (132) 

with the same coefficients b and c as in (130). There- 
fore, the leading finite-size corrections to /5c when de- 
fined by Re - Rk cancel. At criticality, 

Re{/3c) = Rki/ic) = Rc+ cA?;" + ■ ■ ■ , (133) 

so the leading corrections only move the intersection 
point of Re(J3) and Rk(fi) vertically without shifting the 
so defined critical coupling. 

The results for SU(3) from the R^ - Rk intersection 
points are compared to the extrapolated pseudocritical 
couplings from intersecting Rk(/3) with its universal crit- 
ical value Rc via Eq. (131) in Tab. 1. 

In 2-1-1 dimensions the critical couplings grow lin- 
early with A^, to leading order at large A^,, c.f. Eq. (61), 
and the critical temperature in units of the dimensionful 
continuum coupling gj, is given by the slope, 

/3c(N,)l{2Nc) = {Tclgl) N, + --- . (134) 

From the values for f^c with A^, = 4, 6 and 8 in Tab. 1, 
one obtains for SU(3) [100], 

^ = 0.5475(3) (135) 

corresponding to Tc/ yfcr - 0.9938(9) with a zero tem- 
perature string tension y[(rl^^ - 0.5509(4) from a 
weighted average of the four values in [133]. This is 
consistent with T,/ a/o^ = 0.9994(40) from [99]. 

IVloreover, because the spatial center vortex free en- 
ergies Fk for sufficiently large L depend only on L'^''f, 
and t oc (J3 - /5c) at leading order in A^,, when expanding 

Fk{Ji)^- In Rc + diN,)(p -/?,) + ■■■, ( 1 36) 
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N, 


/3c{Rk = Re) 


/3c {Rk = Rc) 


Lit. 


2 


8.15309(11) 


8.15297(57) 


8.1489(31)^ 


4 


14.7262(9) 


14.7194(45) 


14.717(17)"f 


6 


21.357(25) 




21.34(4) t 


8 


27.84(12) 







Table 1: SU(3) critical couplings from selfduaUty (weighted means), 
and intersection with the universal value (extrapolated) from [100], 
previous literature values from ^[99], *[108]. 



one expects the slope at /3c to behave as 

d{N,) ~ . (137) 

Our current best estimates from fitting these slopes are: 

V = 0.99(3) for SU(2) where one expects v = 1 for the 
2d Ising model, v = 0.818(24) for SU(3) where v = 
5/6 « 0.833 for the 2d 3-state Potts model, and v = 
0.673(9) for SU(4) as compared to v = 2/3 for the g = 4 
Potts model. An update on our results including more 
on SU(4) will be pubhshed elsewhere. 

For SU(4), the Z4 center symmetry alone does not 
uniquely specify the effective spin model to describe 
the dynamics of Polyakov loops. SU(4) is a rank-three 
group and has three fundamental representations, 4, 4 
and 6. So even the simplest effective Polyakov-loop 
model will consist of two distinct real terms, with near- 
est neighbor couplings between loops in 4/4 represen- 
tations and between loops in the 6 representation, c.f. 
Ref. [134J. Depending on the relative weight between 
the two, the corresponding spin model could be any of 
the Z4-symmetric Ashkin-TeUer models with three en- 
ergy levels per link and continuously varying critical 
exponents between the q = 4 Potts model class with 

V = 2/3 and yS = 1/12, and that of the planar or vector 
Potts model, or simply the clock model, which corre- 
sponds to two non-interacting Ising models in this case 
with V = 1 and yS = 1/8. In between these limits 
the model corresponds to two interacting Ising models 
with spins s,- and cr, and it thus has two order param- 
eters, the usual magnetic ones or (cr,), which are 
the same in the isotropic model, with critical exponent 
j8m = (2 - >')/(24 - 16^), and an electric one <s,cr,) with 
exponent /3e - (12 - 8^) \ with = for the standard 
Potts model (j3e - /3m), and y = 1 for the non-interacting 
case (fie = 2/im) of the clock model [122]. 

The present conclusion from the studies in Refs. [97- 
99] is that the deconfinement transition in {2+\)d SU(4) 
gauge theory is weakly 1 ^' order. We do observe, how- 
ever, at least approximately at the length scales corre- 
sponding to our spatial lattice volumes, a universal scal- 
ing which seems closest to the standard q = A Potts 



case [100]. The critical couplings from the universal 

amplitude ratios in Eqs. (59) are fully consistent with 
those determined independent of Potts scaling and uni- 
versality, e.g., for AT, = 4 we find pc = 26.294(2) from 
Eqs. (59) as compared to /5c - 26.283(9) independent 
of that. The correlation length exponent is also consis- 
tent with 2/3 and the data scales nicely for all volumes 
up to A', = 80. This maybe still too small to see the 
P' order nature, and the conclusions may well depend 
on A', as well. But even if the transition is weakly 
order in the infinite volume {Ns 00) and continuum 
limits (A^, — > 00), it seems legitimate to ask, why it is 
the standard Potts model universality that is observed 
before those limits are reached, and why not any of the 
other Ashkin-TeUer models, and can we derive the cor- 
rect effective Polyakov-loop model to demonstrate that? 

3.5. Fractional electric charge and quark confinement 

Center symmetry is explicitly broken and the fi- 
nite temperature transition of QCD becomes a smooth 
crossover when dynamical quarks in the fundamental 
representation are included. This crossover at zero 
chemical potentiid is very well studied on the lattice 
[135, 136]. In presence of fields such as the quarks in 
QCD which faithfully represent the center of the gauge 
group, there are no twisted boundary conditions on the 
torus and therefore no center vortex ensembles to study. 
The detailed understanding of the deconfinement phase 
transition in the pure gauge theory might appear to be a 
rather academic exercise then. 

This neglects the electric charge of quarks, however, 
which is commonly expected to only require small per- 
turbative corrections. On the other hand, the inclusion 
of the quarks' fractional electric charge is well known 
to lead to a global Zg symmetry of the full fermion and 
Higgs sector of the Standard Model. This Ze symme- 
try combines the centers of the color and electroweak 
gauge groups, see also [137], and [138] for a review. 
Since quarks carry fractional electric charges Q - |e or 
-\e, their color and electromagnetic phases in the pair 
of combined transformations, 

^^ilnli^^ilnQie^^ ^^-ilnl^ ^^-ilnQle^ ^ g^j^g^ ^ U(l)em, 

cancel precisely. They act trivially on all other particles 
in the Standard Model, which are blind to the center of 
SU(3) and carry integer electric charge. Electric charge 
Q is related to hypercharge Y and the third component 
of weak isospin hy Qje - ti + 7/2, with e'^'' = -1 e 
SU(2) and Y quantized in units of 1/3. The synometry 



L. von Smekal /Nuclear Physics B Proceedings Supplement - preprint (2013) 1^2 



35 



is therefore generated by 

(e'2''/^-l,e'''^) e SU(3) x SU(2) x U(1)k, (138) 

which gives six elements including the identity. 

A global symmetry brings with it the possibility of 
a phase transition characterized by spontaneous sym- 
metry breaking. In this case, one expects the transi- 
tion to be driven by topological defects (vortices) that 
carry both color and electromagnetic flux, that is, cen- 
ter vortices with an additional electromagnetic Dirac 
string [139, 140]. Could the physical realization of this 
symmetry have non-trivial implications for confinement 
and the phase structure of QCD as part of the Standard 
Model then? We have started to address this question 
recently in a toy model with half-integer electrically 
charged quarks in a two-color QCD world with elec- 
tromagnetism [141, 142]. 

The starting point is QC2D plus electromagnetism 
with 2 flavors of Wilson fermions in 3 -1- 1 dimensions. 
By including 'up' and 'down' quarks with fractional 
charges ±e/2 relative to the U(l)em gauge action, we 
obtain a model with a global Z2 symmetry. The lattice 
action is 

5 =-^(^tl-C/^ +y6emCOS0^) + 5/,w, (139) 
□ 

where S f^w is the usual Wilson fermion action with the 
distinction that parallel transporters for quarks are prod- 
ucts of SU(2) color matrices and U(l)em phases, 

^/;,(x)e'■''''W/^ Uf,(x) e SU(2), 0f,(x) e (-27r,27r] . 

(140) 

The SU(2) plaquettes and U(l)eni plaquette angles 
6^ are formed from [/^ and 6^ in the usual way. In this 
model, 'fractional charge' means that the parallel trans- 
porters for quarks contain half the U(l)em angle relative 
to the Ofi's that appear in the plaquette angle 9^. That 
is, an e'^'p^^ = -1 electromagnetic link for quarks ap- 
pears as an e"''' = +1 link in the U(l)em gauge action. 
An important point is that the 'volume' of the com- 
pact U(l)em is determined by the quarks, which carry 
the smallest quantum of electric charge. The range of 
6^{x) is chosen such that we integrate over all possible 
electromagnetic transporters for the quarks, amounting 
to a double counting in the U(l)em gauge action and for 
all integer charged particles. This is consistent with the 
premise that our compact U(l)ein is the result of sym- 
metry breaking in a SU(3) SU(2) x U(l)em unified 
theory, e.g., see Ref. [143]. 

It is clear from Eq. (140) that a color center element 
-1 G SU(2) combined with an electromagnetic phase 
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Figure 24: The SU(2) Polyakov loop in tlie pure gauge theory (green), 
in full two-color QCD with ff = 0.15 (red), and in the SU(2) X U(l)era 
toy model with electromagnetic couplings varying from total U(l)eni 
disorder at ji^^ = to deep into the Coulomb phase at ySem = 2 on a 
16^ X 4 lattice with Wilson quarks. 

g'«',i/2 - _i act as the identity on the quarks. The model 
therefore retains a color and electromagnetic combined 
Zi center symmetry, despite the introduction of dynami- 
cal quarks. This is the analog of the hidden Zg symmetry 
of the Standard Model. 

The results in [141, 142] indicate that the usual order- 
ing of color links by the dynamical quarks is negated by 
the inclusion of fractional electric charge. This is under- 
stood by analogy with the spin systems. In QCD alone, 
terms in a loop expansion of the fermion determinant 
that wind around the temporal direction favor the center 
sector in which the traced Polyakov loop Pcoi = 1- They 
break center symmetry and lead to an ordering in much 
the same way as an external magnetic field H does for 
the Potts model spins in Eq. (50). The inclusion of frac- 
tional electric charge in our model bestows quark loops 
with an additional U(l)en] phase which may undo this 
eff'ect. Consider, for example, the hopping expansion 
of the Wilson fermion action for N, - 4 time slices to 
leading order in the hopping parameter, O(k^), 

(0 
^cos^-trf/^ (141) 
□ 

+ 8 cos iPg,2i^) ■ Re PeoK^O J + ■ ■ • • 

If the U(l)em Polyakov loop angle for quarks Pgjiix) - 
2|^o' 0t{t,x)/2 is disordered, then the Pco\ - 1 sector 
is no longer favored and the SU(2) center symmetry is 
dynamically restored. The efifect is analogous to placing 
a Potts model in a fluctuating magnetic field, or to the 
Peccei-Quinn mechanism [144], in which the coupling 
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of the CP violating term in QCD to an axion field allows 
for the dynamical restoration of CP symmetry. 

Since the parallel transporters for quarks possess a Z2 
degree of freedom that the U(l)em gauge action is blind 
to, this is possible even in the Coulomb phase for integer 
electric charges. Indeed, the SU(2) Polyakov loop in the 
toy model with half-integer charged dynamical Wilson 
quarks is indistinguishable from the quenched result at 
values of the hopping parameter, or the quark mass, that 
would otherwise cause a significant amount of ordering 
in standard 2-color QCD. As shown in Figure 24, this is 
equally true when the U(l)em links are totally disordered 
in the confined phase jSem ^ 1.01, and in the Coulomb 
phase at ySem 1-01 of the U(l)em- Since the U(l)ein 
gauge action is blind to the Z2 c U(l)em disorder as seen 
by the quarks, i.e., e'^f^^ = +1, it is unable to remove it. 

The quenched result in Fig. 24 reflects the second 
order nature of the phase transition of pure SU(2) in 
the universality class of the 3c/ Ising/Z2 gauge system 
with typical finite-size corrections. It clearly turns into 
a crossover in full QC2D at a- - 0.15 as in the Ising 
model with a symmetry-breaking external field. When 
the half-integer electromagnetic charges are included, 
the SU(2) Polyakov loop falls straight back onto the 
quenched curve indicating that there might be true dis- 
order in the vacuum when QCD is embedded into the 
Standard Model. 

This is not the full story, however. A comparison of 
the n and p meson masses in the toy model and in QC2D 
with equal hopping parameter k and SU(2) gauge cou- 
pling yScoi as in [145] reveals that the mass scales have 
changed dramatically in the toy model [142]. One plau- 
sible interpretation is that there is strong Z2 disorder 
even for (5^^ — » 00 which adds to the strength of the con- 
fining potential. This effect might be unphysical. One 
needs to go back and verify scaling of the masses (or 
correlation lengths) in physical units as the continuum 
limit is approached which might not even exist without 
unification in the toy model. 

Meanwhile, we have explored the phase diagram at 
the leading order (9(/c'*) in the hopping expansion for 
A^, - 4, with the effective action in Eq. (141), which 
allows us to cheaply check our intuition about the result- 
ing couplings of SU(2) and U(l)em links, without hav- 
ing to worry about the chiral limit for k. We thus take 
the leading-order hopping expansion at face value, as a 
model beyond its range of applicability as an approxi- 
mation to the full toy model. This effective model then 
shares its qualitative features with a fractionally charged 
fundamental Higgs model [146]. 

Fig. 25 shows the results for the SU(2) Polyakov loop 
in the effective SU(2) x U(l) model with the leading or- 




Figure 25: The SU(2) Polyakov loop on the 16^ X 4 lattice with the 
leading order interactions from Eq. (141) in the /3e„, — ► 00 limit from 
Ref. [142]. 



der interactions in Eq. (141), where the U(l)em phases 
g'fl^(-f)/2 ^j-g restricted to +1, which amounts to Pem °° 
in the U(l)ein gauge action. As k increases, the disorder- 
order deconfinement transition of the SU(2) Polyakov 
loop moves to smaller values of yScoi and sharpens dra- 
matically. Note that in the combined limit k, Pem — > °°, 
the plaquette-plaquette coupling in Eq. (141) forces the 
SU(2) plaquettes to take the values +1. For very large 
K the transition line should therefore terminate with the 
first order bulk transition of the Ad Z2 gauge theory at 
yScoi ~ 0.44 [147]. As the SU(2) plaquette is driven to 
unity for large jScoi, the U(l)em plaquettes correspond- 
ing to quark loops, cqs{9^ /2), receive an effective cou- 
pling that suppresses Z2 disorder. The relevant U(l)em 
Polyakov loop cos f e/2 remains disordered at small val- 
ues of K, but transitions to unity for large p^oX and k (not 
shown). Due to the Z2 disorder transition, the model 
might have three distinct phases, the confined phase for 
small K and small p^ou the deconfined one for small k 
but Pc < Pco\ and a Higgs phase in the corner where both 
couplings K and p^ox are large. We currently investigate 
the same effect of Z2 disorder from half-integer elec- 
tric charge in the phase diagram of fundamental SU(2) 
Higgs models [146]. 

Given the fractional electric charge of quarks and the 
existence of a global symmetry that relates the centers 
of the color and electroweak gauge groups, it may be 
misleading to study non-perturbative phenomena such 
as confinement and perhaps also dynamical mass gen- 
eration in QCD alone. In the toy two-color model, 
the coupling of Wilson quarks with half-integer elec- 
tric charge relative to a compact U(l)em has a dramatic 
effect on the color sector This is not due to electromag- 
netic fluctuations which should have small effects deep 
in the Coulomb phase for large ySem. but it can be due 
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to a frozen Z2 disorder that the low energy theory inher- 
its from the grand unified theory. Such frozen disorder 
could have an effect akin to a random external field in 
Potts models [148] which are known to have phase tran- 
sitions where they would not otherwise have them, with 
constant extemal fields. 

Here the SU(2) Polyakov loop shows an order- 
disorder transition that is consistent with the sponta- 
neous center symmetry breaking transition of the pure 
gauge theory at values of the hopping parameter where 
center symmetry is clearly explicitly broken in stan- 
dard two-color QCD. This happens despite the fact that 
static quarks are now described by the combined SU(2) 
and half-integer U(l)em loops which continue to show 
a crossover behavior. The scales of both transitions 
are still widely separated and this problem needs to be 
solved in such models. 

Combined twisted boundary conditions are neverthe- 
less possible with more than one gauge group and the 
corresponding center vortex ensembles might become 
relevant for a deconfinement phase transition in the stan- 
dard model if its hidden center-like symmetry is rele- 
vant to the real world. Whether this is the case or not 
seems hard to decide without grand unification. Mean- 
while one might want to find out whether SU(3) x Z3 
models of the kind of our two-color toy model in the 
ySem ^ °° limit are at all suited as effective low energy 
theories for hadron physics. If we did have absolute 
confinement in the real world, separated by a true phase 
transition from the quark-gluon plasma, this would ob- 
viously make a conceptually and phenomenologically 
important difference as compared to our present under- 
standing that quarks are only exponentially suppressed 
in a hadronic phase that is analytically cormected to the 
high-temperature plasma. 

4. Concluding remarks 

Many interesting facets of strongly interacting mat- 
ter can be studied in QCD-like theories. Changing 
the number of quark flavors, their masses or the num- 
ber of colors often leads to simphfications or idealiza- 
tions which allow to understand qualitative features and 
mechanisms from such powerful concepts as universal- 
ity, scaUng and finite- size scaling. They allow to re- 
late phenomena from diverse areas of physics by often 
surprising analogies. The hquid-gas transition of nu- 
clear matter, Heisenberg ferromagnets for chiral sym- 
metry breaking, ultracold fermionic quantum gases with 
BEC-BCS crossover for two-color QCD at finite baryon 
density, the very well studied spin models and their du- 
ahties for the deconfinement transition and the electric 



fluxes in the pure gauge theories are some of the famous 
examples touched upon in these lecture notes. 

The description of critical phenomena in statistical 
physics and quantum field theory is today synonymous 
with the renormalization group. Valuable introductory 
accounts of various aspects of the renormalization group 
have been given at this 49* Schladming Winter School. 
I have presented some selected examples of applications 
of basic renormalization group concepts towards a bet- 
ter understanding of chiral dynamics at finite density, 
quark confinement, the different phases of strongly in- 
teracting matter at extremes of temperature or density, 
and the nature of the transitions between them. 

Arguably, two of the most promising non-perturba- 
tive tools to study the QCD phase diagram at the mo- 
ment are the functional renormalization group and lat- 
tice simulations, especially in combination. These two 
complement each other very well. Non-perturbative 
functional methods in general rely on additional as- 
sumptions. These can be guided by intuition, derive 
from fundamental requirements such as causality, local- 
ity or gauge invariance, or they can be backed up by 
results from other approaches as external input. This 
is what lattice simulations can provide where they are 
possible without fermion sign problem. Therefore, it 
will be extremely valuable to investigate and ultimately 
completely understand the phase diagrams of QCD-like 
theories such as two-color QCD or the G2 gauge theory 
which don't have this problem. 

A lot of work in this direction has been done in two- 
color QCD as reviewed in Section 2. Here we are natu- 
raUy interested in the functional renormalization group 
treatment which we have seen to be well capable of 
describing the effects of finite baryon density and the 
competing dynamics of collective mesonic and baryonic 
fluctuations. The bosonic nature of the baryons in two- 
color QCD with a BEC-BCS crossover provides inter- 
esting analogies, but I have argued that the impact of 
the finite baryon density on the chiral transition reveals 
a genuine effect. There are reasons to believe that there 
might not be enough chiral symmetry breaking left for 
a chiral phase transition and a critical endpoint at fi- 
nite baryon density in the QCD phase diagram either. 
Suitable theories with fermionic baryons but without 
fermion sign problem seem to be the logical next step 
to study at finite density. The G2 gauge theory is such 
a theory and the studies of its phase diagram have only 
just begun. The challenge there will be to disentangle 
finite density effects due to bosonic from those due to 
the fermionic baryons, an ideal task for the next com- 
bined effort with the functional renormalization group 
and lattice Monte-Carlo simulations. 
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The confinement problem is obviously a hard one. 

But a better understanding on a fundamental level even 
in the simple 2+1 dimensional SU(A') gauge theories 
from exact results seems worthwhile to work our way 
up from there. I found some of these results such as 
the selfduality reflected in the center vortex and elec- 
tric flux ensembles quite enlightening and perhaps also 
unexpected. The models may not be realistic but the 
results are quite powerful, and who knows, perhaps we 
will eventually find that it wasn't wasted but that com- 
bined vortices and disorder in the vacuum may be rele- 
vant when QCD is embedded in the Standard Model. 
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